'Digital simulation and optimization 

OF CATALYTIC REFORMING UNITS 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


By 

ASHOK KUMAR SHARMA 


to the 

DEPARTMENT OF CHEMICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

MAT 1973 




l.l.T. KANPUR 
CENTRAJ- yBRA^ 

^ NO. 



fiJ MAR 


(^Q>C) ^ Si 




He - ni- 3"'^'- SHW" 



il 


CERTIFICATE 


This is to certify that the present •work 
'Digital Simulation and Optimization of Catalytic- 
Reforming Units' has been carried out under ray super- 
vision and that this ■work has not been submitted 
elseiAtoere for a degree . 


May 8, 1973 





(Ml Gopala Rao) 
Associate Professor 
Department of Chemical Engg. 
Indian Institute of Technology 
Kanpur-208016, India 



Aomo v^SDSmiBlT s 


The author wishes to express his deep sense of gratitude 
to Dr. M, Gopala Rao, for his ■valuable guidance, conunents, 
criticism and encouragements throughout the development of 
this thesis. The author is also highly indebted to Dr.A.K. 
Wagle for his help and keen interest in developing the 
optimization routine . 

The necessary data provided by Engineers India Limited 
and Indian Petrochemicals Corporation Ltd. , is graijeifully 
acknowledged. 

Thanks are due to Mr. J.K. Misra for typing this thesis. 

I* 

Research assistantship provided. 'by ChemJ-cal Engineering Dept, 
is highly acknowledged. 


A. K. Sharma 



ABSTRACT 



A fiiathematical model for the. simulation of catalytic- 
refoming units has "been developed. The model takes into 
account most of the complex reactions oc curing in the process 
The model vra.s used to simulate on a digital computer two 
industrial .Platforming units. Agreement between simulation 
results and plant data is very satisfactory. It has been 
found that the performance of an industrial unit can he 
improved significantly hy changing some of the present operat 
ing conditions. An optimization routine to minimize the 
weight of catalyst and the operating cost has been developed* 
Piuther work can he done to combine mathematical models of 
inter heater and a flash unit with the present simulation 
model. 
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MSTMCT 


A imatbematical model for the- simulation of catalytic- 
reforming units has been developed. The model takes into 
aceount most of the complex reactions- occuring in the process 
The model was used to simulate on a digital computer two 
industrial .J?latfoiming units. Agreement between simulation 
results and plant data is very satisfactory. It has been 
found that the performance of an industrial unit ean be 
improved significantly by changing some of the present operat 
ing conditions. An optimization routine to minimize the 
weight of catalyst and the operating cost has been developed* 
further work can be done to combine mathematical models of 
inter heater and a flash unit with the present simulation 
model. 
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CHIPTER I 


DJTRODUCTIOl 

Use of analogy and simulation has contributed a lot 
to progress of science and technology. Through analogy with 
familiar and simple events one can understand more complex 
phenomena of the world in which one lives. To test analogi- 
cal conclusions, experiments have to be conducted which 
simulate the natural conditions one wishes to understand. 
Because c£ its very definition, an analogy or simulation can 
not be identical with realities of nature it attempts to 
represent . 

As technology expands, the demand for laboratory 
simulation also expands. In some experiments, the parameters 
which govern the behaviour of the system under study are 
complex enough to restrain the experimenter from such a work. 

■As a result modern research and technology have found digital 
and analog simulation not only expedient but also indispens- 
able. And now with the advent of big and fast computers, one 
can manipulate the mathematical equations of the system, called 
mathematical model, with a reasonable degree of effort and 
cost. 

Digital simulation is also of vital importance for the 
chemical industry as it can be used; 
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1. To facililjate design and operaiiion of process or 
its parts. 

2. To choose operating- conditions for changed production 
pattern or increased production in an existing plant. 

3. To study the effect of extreme ranges of operating 
conditions, some of which might be impractical or 
impossible to use in a pilot plant. 

4. To compress or expand real time, 

5. io compare various proposed designs and processes 
not yet in operation. 

6. To test the hypothesis about a system or process 
before making a decision. 

7. To study the effect of changes in variables and 
parameters with reproducible results. 

8. To gain sufficient knowledge about the system to 
facilitate computer control. 

It might appear from the discussion presented above 
that it is much better to analyse the system by digital- 
simulation than to analyse it by carrying out laboratory 
experiments or by pilot-plant studies. Unfortunately, there 
are certain limitations to the simulation techniaue as the 
accuracy of the results obtained by it depends largely ;on 
correctness of the mathematical model developed and the 
accuracy of the physical and chemical data that go into the 
model. 
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In India "to meet the increasing demand of aroaiatics 
(installed capacity of benzene in India has incr^asad fro® 

15,000 tons per annum in 1963 to 87,000 in 1969/70 and the 
requirements for 1973/74 have been estimated at 1,40,000 tons) 
a catalytic-reforming plant is being commissioned in Gujrat 
State. At present the process know-how for this process whfcla, 
uses naphtha as feed stock, is imported from the Western 
European countries. Catalytic-reforming when combined with 
solvent extraction or extractive distillation produces- , 
benzene for ultimate conversion to phenol or styrene, tolus-Xie 
for nitration purpose, xylene for plastics etc., in qijanti- . 
ties far beyond the capacity of any other source to provide 
these materials. Catalytic - reforming whioh goesi, 
various names such as Platforming (Universal Oil Product Co.) , 
Power- forming (Esso Research and Engineering Co,), Ultra- 
forming (Oil Co. of Indiana), etc. is also used to upgrade 
the octane number of straight run gasolines, naphtha and 
other feed stocks. The hydrogen which is liberated in this 
process can be ploughed back into heavier portion of crude 
oil either by direct hydrocracking of heavy oil or by milder 
hydrogenation to produce superior Diesel fuel or improved 
charge stock to the catalytic cracking units. In order to 
achieve self-sufficiency (partially or fully) in this important 
process, it has been intended in the present work to simulate 
the reactors system coming up in Gujrat State using the IBE/l 
7044 computer. 
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Calialytic reforming is a bigh pressure arid high ieinpe-- 
raiure process wnich uses balogenaied alumina coaied wiih 
platinum as catalyst. Generally, three to five steel reactors 
packed with catalyst with an inter-heater before each are used 
in series. The reactor vessels are insulated by monolithic 
ceramic liners, light gases which have about 85 percent 
hydrogen are seperated in flash unit from catalytic refonaate. 
A part of this gas is used as recycle gas. Detailed process 
flow diagram is given in kig. 1. 

In the present study, a mathematical model has been 
developed to predict the product stream of a catalytic- 
reforming reactor system for a given feed composition, 
catalyst distribution, inlet temperature, inlet pressure and 
recycle ratio. The effects of temperature, pressure and 
recycle ratio have been studied. This mathematical model 
together with principles of calculus of variation, are used 
to write an algorithm to extremize an objective function 
with constraints on final product stream. Using this algori- 
thm, respective objective functions involving total weight 
of catalyst, and operating costs have been minimized for a 
three, four, and five reactors set-up. 

Present study deals only with the steady state 
behaviour of the process, dynamic response has not been 
studied . The composition of recycle gas from the flash unit 
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is not oaloulated, and instead a fixed composition is 
assumed. The effect of temperature^ pressure and impurities 
in naphtha on the catalyst life is not considered. 




CHAPTER II 


LITERATURE OH GATAIYTIG BEFOMIHG- REACTION 

]\iany cbemical reactions, some simple and some complex, 
constitute catalytic-refoimiing but the reactions leading to 
the formation of aromatics from naphthenes and paraffins are 
of greater importance. Under the catalytic-reforming condi- 
tions, dehydrogenation of naphthenic hydrocarbons is very 
rapid and contributes substantially for the yield of aromatics. 
These reactions are highly endothermic and result in somewhat 
to volumetric shrinkage since the aromatics product is denser 
than the parent naphthenes (1). Aromatics are also foamed 
from naphthenic intermediate resulting from cyclization of 
paraffins. But this reaction is not obtained in most ideal 
sequences because even at a moderate severity level, naphthene 
dehydrogenation is close to completion while paraffin dehy- 
drocydization is still in the initial stages, furthermore 
for all types of catalyst used commercially, the dehydrocycli- 
zation of paraffins is less efficient than dehydrogenation 
of naphthenes, for this reason, naphthene dehydrogenation 
to aromatics is considered to be the backbone of catalytic- 
reforming (2). Hydrocracking which is opposite of dehydro- 
genation from the stand point of yield and the heat of reaction 
is also an important reaction. In this reaction high-molecular 
weight hydrocarbons decompose to lower-molecular-weight 
hydrocarbons. 
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As discussed above, catalybic-refonaing is simultane- 
ous dehydrogenation, cyclization and hydrocracking* One can 
imagine the complexity and cumbersomeness if rate expressions 
written for each carbon number compounds in the naphtha 
undergoing the reforming reactions, are to be solved* Kinetic 
studies have been done by taking each hydrocarbon seperately. 
Attempts have also been made to simplify the problem by treat- 
ing naphtha to be of some average carbon number depending 
on its composition. In such studies Smith (3) has come cut 
with a reaction model and rate equations #iich are reported 
to be applicable for naphtha feed stocks under widely used 
operating conditions for catalytic-reforming. Katrush and 
others (4) report that these kinetic equations fail to 
predict plant behaviour for Russian naphthas. Ror each 
reaction, the rate equations of each reaction used by Katrush 
et al, for three reactors in series have different rate cons- 
tants for each reactor.. Their results are based on the 
hypothesis that in the first reactor the main reaction is 
dehydrogenation? in the second dehydrogenation and cyclization, 
and in third dehydro cyclization, isomerization and hydro- 
cracking. Recently Heningsen and lielson (5) have came out 
with different reaction model and kinetic rate expressions. 

But these have been reported only for Cg system and reactions 
occuclng at 500°C and 30 atm. The rate expression of Marianme 
Vermorner and others (6) who carried out their work seperately 
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are approximately the same as obtained by Smith. Smith ( 3 ) 
bas proposed the following reaction model: 

1* Naphthenes ^ /romatics + Hydrogen 

2. Naphthenes + Hydrogen ;;frr Paraffins 

Cv, + H„ r w 

^ 2n ^2 ^ 2^+2 ( 2 , 2 ) 

5. Paraffins + Hydrogen -* Hydrocarbons 

^n ®2n+2 + “2^ ^2 J 5 + ^ 2 ^ 6 '^ •**%^12^ ^2.3) 

4. Naphthenes + Hydrogen Hydrocarbons 

(CH^ -f CgHg^ (2.4) 

Ir the above reaction model, it has been assumed that 
the hydrocracking of aromatics is too slow and can be neglected. 
The stoichiometry of reaction (2. 3 )and( 2. 4) which shows that 
equal number of moles of methane ,, ethane , propane,' butane and 
pentane are formed from hydrocracking of naphthenes and 
Paraffins, has been derived from experimental observations, 
that these species are formed approximately in equal molar 
proportion. 

Prom thermodynamical considerations, one expects 
catalytic reforming reactions to proceed at low pressures 
and high temperatures. But at high temperatures, coke gets 
deposited on the catalyst due to hydrocracking, v\hich reduces 
the catalyst activity. Coke deposition can be overcome by 
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having a high hydrogen partial pressure which can he obtained 
at higher operating pressures. It also reduces sulfur deposi- 
tion on the catalyst by converting sulfur to hydrogen sulfide (7) 
Operation at high pressures lengthens the tine between regene- 
rations, increases 'yield of aronatics and reduces the amount 
of light end products by suppressing hydrocracking. But as the 
pressure is increased beyond a certain limit, the yield of 
aromatics deteriorates because of more hydrogenation of naphthe- 
nes to paraffins. 

Hydrogen partial pressure can also be increased by increas 
ing the hydrogen to hydrocarbon ratio which is direct function 
of the recycle ratio. Though at low recycle ratio more hydro- 
gen can be formed, the partial pressure of hydrogen is not 
sufficient to suppress the coke deposition on the catalyst. 
Increase of hydi'Ogen concentration can also reduce the yield 
of aronatics because it reverses the order of desirable reac- 


tions. 



CHilPTER III 


MiiTHEMAIIGAL MOI ELl i& AND SBtULaTIOE 


Asstming that hydrocracking of aromatics to hydro- 
carbons is negligible and methane, ethane, propane, butane 
and pentane are for.ied in equal molar proportions from the 
hydrocracking of naphthenes and paraffins, one can write 
the over all reaction model for catalytic - reforming reactions! 


1 . 


2 . 


3. 


4. 


Naphthenes ^ Aromatics + Hydrogen 

3Ho 






Naphthenes + Hydrogen Paraffins 




Paraffins + Hydrogen 


Hydrocarbons 


“ 211+2 ¥ “2 


i (CH^ + OgHg + . 


(3.1) 


(5.2) 


(5.3) 


Naphthenes + Hydrogen Hydrocarbons 
°n “ 2 n ^“4 +? 2 V 


(3.4) 

Purthermore, it has been found reasonable to assume that 
each of the three hydrocarbon classes eg. paraffins, naphthenes 
and aromatics in naphtha, has on the whole the sane nimber of 
carbon atoms ’n*. Based on this, value of n can be obtained, 
by expressing molecular weight of naphtha in terms of molecular 
formula as shown below: 
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^ ^2n+2^ + XN* (C^ + Xi. (C^ ^2n-6^ ~ 

03: n = -^ (MWj, - 2 XP + 6 X.i) 

'Experimental data shows that heats of reactions are 
nearly independent of component molecular weight over the 
range normally encountered in reforming process, when hased 
on moles of Hg entering into the reaction. Based on this, 
kinetic rate expressions and heat of reaction for reactions 
(2. 1) through (2. 4) as reported hy Smith (3) are given helowi 

1* = k^ ^ ) 


2 . 


k 


K 


1 


1 

AH 


1 


k^ 


Kg 


AH. 


Ih mole of nanhthenes converted to aromatin-q 
tr (lb of catalyst)"' 

C3.5) 

exp (23.21 - 34750A) ^ 

exp (46.15 - 46045/1) atm^ 


= 30500 


Btu 


16 mole Hg produced 


(Bw B 


P. 


P 


2 “ K. 


) 


lb mole of naphthenes converted to naraffins 
~~ (hr. lb cat.) 


exp (35.98 - 59600/1) . ai-gg . ls .. nwhthenes 

hr. “ 

exp (8OOO/T - 7.12) atm" 


hr. lb cat atm*^ 
„-1 


Btu 


■ 19000 mole Hg consumed 


0<6) 
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3. 


4. 


k 


3 

AH. 


"■4 

AH. 


= k 


■-z -V— 

J -tfJJ 


C5-7) 


= exp (42.97 - 62300/T) 


k 


24300 — ^ 

mole H 2 consumed 

•D 

'H 


4 


(3.8) 

exp (42.97-62300/T) ^hr° ' ° ( ^ 


nole Hg consumed 


Under catalytic- reforming conditions for T/vbicla tempe- 
rature and pressure are of the order of SOO-IOOO^P and 15*»50 
atmospheres respectively, the reaction system is always in 
gaseous state. Preliminary calculations showed that the 
flow of these gases in the reactors is in highly turbulent 
region with of the order of 15,000 or more. Since there 

is sufficient packing on both the sides of reactors, it can 
be assumed that the flow of gases is in fully developed tur- 
bulent flow region in the catalyst zone. Catalytic-reform- 
ing reactors also have L/d'.p » 20 and d-^/dp > 10. Por such 
flow conditions and reactors one can safely assume a plug-flow 
pattern (S) . 


If ^i3 ^14 are extents’ of reactions (2. 1), 

(2,^, (2.3) and (2.4); per mole of total feed at a place in 
reactor i, then the partial pressiire terms appearing in the 
rate equations can be written as follows: 
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P. 

ij. 


Piji mole fraction of aromatics 
Pr, 

1 'xorax mox6 

XI.I + X. 


= P 


T 


== Pi 


T 


•K 


•f 


■f 


moles of aromatics 
Total moles 


JJ. 


1+3Xi^-Xi2+(i^.5-1- ^ )X^j+(“^.5-1- 


X.1I + X . , 
11 


1+3X,^-X.2-X,^ 


XSI - (X,,+X,,.^X,,) 
1+3X.,-X.2-X,^ 

1+3X,^ - ^12 - Xl4 


XHI + 3X. ^-X.2- Y ^13 - T 


i4 


'H 


•T 


1+3X,^ - X,2 - X,^ 


Material balance for reaction fe. 1) for a differential 
reactor section of catalyst weight awj gives; 


or 


hXi^ ' 

dX. . 

i1 

dw 


■T 


r^ AW 


T 


T 


( P 


T 


(XHI + x^^ + x.,^)) 




Xi2 - X,^ 


K^(l+3Xj_^ - Xj_2 - Xj_^)^ 


T 

(3.9) 
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For reactions (2.2), (2.3) and (2.4) also similar 

relations are obtained? 

,2. 


dX 


i2 


dw 


^2^ 


Pj(XNI-X.^-X.^^X.^)(XHI-^3X.^-X,,jy fc.^) 

(1+3X^,-X.2-X.^)2 


P^(XPl4-X.g-X.^) 

■ K2(1+3X,^.X,2-X,^) 


) / P 


T 


( 3 , 10 ) 


dX 




dw 


5 (1h-3X.^-X.2-X.^).Fj 


(3.11) 


dX. , 
i4 

dw 




(xm-x,,-x,,-x,^) 

(1.3X,^-X,2-X.^)P, 


< 5 . 12 ) 


Since reactions are carried out in an adiabatic 
reactor system, energy balance equation can be written as: 

Cp.Fij.AT = ( 4 X^^. 3 (-AH^) + 4 X. 2 (- AH^) 


+ 4X^^ ^ (-AH^)+AX^^ I (-AH^)).?^, 


Dividing both sides of the above equation by aw, 
rearranging, and taking the limit Aw - 0, one gets, 


dT 

"Kw 


dX. . dX_ 

( 3 (-AH^) + -.rr^ (-AH 2 ) + 


^ a|2 (..Hj) + ^ I (-AH^)) / Oj, (3. 13) 
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. After manipulating the ec^uation normally used for 

I 

packed beds (9), pressure drop in the differential element 
can he written as follows: 

® • 75 ) (1-e) 3X, p 

^ \ 14 ; 

p^.Area^.d 


Initial and boundary conditions for reactor system are: 

0 j = 1,2, 3, 4 (5.15) 

3 = 1 , 2 , 5, 4 ( 5 . 1 «) 

i = 2 , 3 , ♦ . • 

Equation (3.9) through (3.14) together with (3*15) and 
(3*16) constitute the mathematical model for the reactor sys- 
them* 

SHAULATIOH: 


Since fresh feed is mixed with the recycle gas before 
it is fed to the reactors, the solution of equations (3.9) 
through (3.14) can not be obtained without knowing the recycle 
gas composition. Calculations for recycle gas composition are 
difficult because of two reasons: 

Firstly adequate equilibrium data to calculate the 
composition of the gas separated from product stream of the 
last reactor in the flash unit are not available. 

Secondly some iterative scheme has to be used because 


of the recycle stream. 
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The prohlem can he simplified hy asstmiing recycle 
gas composition , "to he the one, noinally encountered in 
catalytic reforming plants. This is justified because of 
the following reasons: 

1. Recycle gas has about 85 percent hydrogen on 
molar basis. 

2. Major portion of the hydrocarbons in the recycle 
gas is constituted by and lower hydrocarbons, 
which do not take part in the reaction model 
assumed for the present study. 

Preliminary calculations also show that a small change 
in tVie recycle gas composition does not cause any appreciable 
effect on the product stream from the reactors. 

Specific heat, C , occuring in eq.uation (3.13) cannot 
be treated as constant because the reactions are carried out 
under adiabatic conditions, and temperature changes are 
large. Specific heats can be assumed to be linear functions 
of temperature and it has been found from preliminary calcu- 
lations that any improvement in this regard does not improve 
results of simulation to any appreciable extent. 

In the specific heat calculations, it was assumed that 
the carbon-number component break-up of naphthenes, paraffins 
(higher than Gg) and aromatics in the reaction mixture is the 
same as it is in the feed to the first reactor. This 
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assumption was made necessary due to lack of information on 
the composition (in terms of carhon numbers) of naphthenes, 
paraffins and aromatics in the product streams. For simpli- 
city, they are lumped together as three species, paraffins, 
naphthenes and aromatics. 

If specific heats of aromatics, na phthenes , paraffins, 
hydrogen, hydrocarbons coming in feed and hydrocarbons pro- 
duced by hydrocracking can be represented by the following 
linear functions of temperature, (details of vAiich appear 
in Appendix C); 

%A " +FA(2).T 

Sn " +FF(2).T 

V = +FP(2).T 

C ,, = FH(1) + FH(2) .T 

pJti 

%E0 = FHC(1)+ FHC(2).I 

C = FHCI(1)+ FIICI(2).0: ; 

pHCI 

the specific heat, C , of reaction mixture can be written as, 

ir 

°p = 0p,_+(MI-X.^-X,2-X.p Cpj, 

+ (XPM.2-Xj_3) Cj,p 

+ (XHI+3X^^-»Xj^2 m ^i3 “ 3 ^i4^ ^pH 

+ (| + I ^14) %HCI 
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2 

^ [XAI,PA(;j)+XPl.J'P(3)+2NI.I'N(j) 

D=1 

+ XHI.PH(;i)+XHCI.FHOl(j) 

+ (Pi(3)-SII(3)+3.FH(j)) 

+ (J5(3)-M(j)-M(3)) X^2 

+ (rHOCj). I - FP(3).PH(3). B^) X^j 

+ (fHC(3) I - M(3) - PH{o) |) Xj^^ ](T)3-'' 

Pressure drops along tbe reactors and inter-beaters 
in catalytic-reforming plants are of tbe order of .7 to -1.7 
atmospheres. These are small in comparison to the operating 
pressures. Preliminary calculations show that the results 
of simulation obtained by assuminis? the nressure as invariant 
are approximately the same as those obtained by considering 
the pressure-drop along length of the reactor. Based on this, 
equation (5.14) was dropped from the system model i.e. dp/dw 
has been put equal to zero. 

Using the recycle gas composi'tion given in Appendix A, 
the system equations (5.9) through (5.13) together with (3.15) 
and (3.16) were solved on the IBM. 7044 digital compnter by 
using the fourth-order Runge-Kutta method. 
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OPTBOZATIOT 

OBJECTIVE : 

Since througliput for a catalytic reforming plant 
is generally very large, even a sna,ll shift in the operating 
conditions and various parameters of the reactor system from 
optimum values can change the economic pattern of the process. 
Reactor sizes, weight of the catalyst, heat exchanger duties 
and compressor size constitute major share of fixed capital 
investment^ while catalyst life, heat duty and power in-puts 
to the compressor, determine major part of the operating 
costs. The variables, temperatures of inlet streams, recycle 
ratio, operating pressure and number of reactors determine 
both fixed and operating costs. Higher temperature of the 
inlet stream though giving rise to higher catalytic activity 
initially, increases hydrocracking lAtoich in turn decreases 
catalyst activity because of the coke deposition on the 
catalyst. This can be overcome by having high hydrogen 
partial pressure vdiich is determined by the choice of the 
operating pressure and the recycle ratio. 

The choice of the operating variables can be subjected 
to various objectives such as minimization of operating cost 
or weight of the catalyst for a fixed throughput and fixed 
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quality of product, or for maxinization of -venture profit 
vtoen installing a new plant. Generally constraints on fhe 
product are there to ensure certain production rate of 
aromatics and a certain concentration of these in the 
reformate. There are also constraints on the lowest and 
the heighest -values of temperature which are determined 
hy reaction rate considerations and catalyst life. 

Mathematically, e,the objective function can he func- 
tion of catalyst weight, catalyst distribution in different 
reactors, catalyst life, furnace sizes and duties and com- 
pressor size and duty. 

The constraints are the following: 

glCXjj^, %3> ^ 4 ) = 

% 2 ’ ~ ^2 
1260 < T^ < 1460 

In the present study, minimization of catalyst weight 
and minimization of operating cost are the two objective 
f-unctions seperately investigated. Catalyst life is a 
variable, which is difficult to take into account. Por 
example the operational life of catalyst depends on tempe- 
rature, pressure, recycle ratio, and feed composition. Very 
little data are available in literature. In the present 
study it is assumed that catalyst is replaced after every 
10.8 months. Also, effect of pressure and recycle ratio 
-was not considered in optimization routine. 
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WiTHmUTICAL lEmOPimT : 


The system under study has reactors in series. Each 
reactor is represented "by five differential eq,uations and 
is related to the reactor next in the series by four equa- 
tions. Extremization of an objective function e f for such 
a system 'which has constraints also on the final product 
stream can be done by usina Calculus of Variant ion and Pontry— 
agin principle o f optimality (10,11,12,15,14,15). 

If is weight of catalyst in reactor i, the reactor 
system can be described by the following equations: 


or 





( 4 . 1 ) 


corresponds to extent of reaction(2,.j)in reactor 
i for 3 = 1 ,2 , 5 and 

X^^ corresponds to temperature in reactor i. 

The initial and boundary conditions for (4.l) are, 
X|^(0) = 0 j = 1,2,...4 , (4.2) 

since reactors are in series 




^i+1,j^°^ 3 = 1,2, 3, 4 (4.3) 

i =1,2,... E-1 


Because of small changes in decision variables the variational 
equat ions of the system can be constructed as, 
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5 9fi- 

1 0 'y ^ + tiigher order terns' (4.4) 

k =1 ik 


Similarly, variation in the objective fxinction , which 
nay depend on the weight of catalyst and various other 
parameters of the inlet and outlet stream at each reactor 
can be written as, 


6 ? 


N 

r’ 

= I 
1=1 


+ f 
1=1 


k=l 

a e 

cW^ 


^ dX^(W^) ^^lk^%^^ 

6y/j + higher order terns (4.5) 


Similarly, variations in the values of constraints, 
which depend on the outlet stream from the last reactor are, 

Sg. = I + higher order terns (4.6) 

^ k=1 ^^Nk 

v" 

6go = Z fiY — + higher order terns (4.7) 

^ k=1 ^Nk 

In the further development, higher order terns have 
been neglected for the sake of simplicity. 

It is readily seen that equations (4.2) and (4.5) for 
( 4 . 4 ) become, 

6X^^(0) =0 3~^j255,4 (4.8) 

, — ,N-1 




(4.9) 
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I'Tow nultiply Equation (4.4) by arbitrary continuous 

functions .(w) and integrate for w = 0 to w = W. +6W.* 

11 

Thus 5 


W, +6¥, 

1 1 

I 

0 




W^+6W- 

/ 

0 


k=l 


«ik 




dw (4.10) 


Integrating left- hand side of eqn. (4.10) by parts. 


W^+6W^ 


I “ / X.. 6X. . dw 

0 13 ID 0 0 ^3 ID 


X, . 6X. . 
ID ^D 


w. +& w . 

1 1 


0 


X,. 6X..(w,46W.) 


ID 


(4.11) 


X,.(0) 6X,.(0) 


, dX,._.CW,) 


ex..(w.) 

[ 61,. (wp + 1 


q^(o) 6X^. (0) 


hj (Wi) eXj^Uwp + 


ir 1' 


6W. + 4^(q)5W.5X.gffp 


13 


13 


H. >,.(wp f. .{wp iW. sw. 


- \3(°) “ 13 ( 0 ) 


(4.12) 
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Neglecting second order terms in (4.12) and making 
use of (4.11), one gets from equation (4»10) 


1 . 

2 . 


Wi+6V/i ^ 

J Y —id 


0 


Wi+6W^ 

I ^^ik ^ii 

k=l ‘^^ik =lk; 10 Q 


^3 SXi- dw 


(4.13) 


low add following to equation (4.5): 

Equation (4.13) for i = 1 , 2 ,..N and 3 = 1 , 2, ...5 

Multiplications of right hand side of equations (4.6) 
and ( 4 . 7 ) hy arbitrary constants and 

3 . Multiplications of equation (4.9) with v. . for • 
i = 1 , 2,..., H -1 and 3 = 1 , 2 , 3 , 4 . 

If and ^2 fixed values and C 2 respectively 

i.®. , 6g^ = Sgg 

= il/ ^ ^ hkOfi)) 


6e 


I~1 4 5e 

* 1J1 kli ^ ^ 


(T.'P 
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5 

I 

k=l 






ae 


"N 
F-1 


®g2 • 


^ ^%k^%) - + yv--71i; \ )• 


N-1 4 

- I J, (^Ik^O) ^ 

1=2 k=1 


Mi5(Wj_) 
0e 


'"i-i.k ■" ax 


Ik 


). 


N 

6X,k(0) - (X^^CO) - ^3^ ) 6X,5 (o) 


4 


^15' 


^ “ ex ., ( O) ^ 


T*' 'Si ^ # 


1k 

5 9f 


Z I / (^iv+ I “Hx ) (4.14) 

=1 k=1 0 ra=l '^^Ik 


There are certain terms which are not at one’s dis- 
posal, namely, '5X^^(’.?j^) in the integrand, 5X^^(0), i = 2,..11, 
3 = 1,..4 and 5Xj[^(T4), i =1,2, 3, 4, 3 = 1,2,, ,5. Therefore 
eliminate these terms from the expression for 6e , by remov- 
ing some of the arbitrariness and making them to satisfy 
the following eq.Ucations: 



i 

k=1 


af 


^ik 


ik 


ex 


= 0 


la 


1.2.. .H 

1.2.. .5 


(4.15) 




g e 


= ° 


i = 1,2,. .11-1 
3 = 1,2, 3, 4 


(4. 16) 



27 


= 0 i = 


2 — 15*#*4 (4*17) 


= 0 

i5 

^N3(%) + ■'l 81§^ + ->2 


i = 1,2,. ,1-1 (4-^18) 


9g. 


Qc 


d = 1;2,..5 (4.19) 

ComlDining equations (4.16) aud (4.17) 


'«‘V ■ 


i = 1,2 ,. .. 1-1 
d = 1,2, 5, 4 


(4.20) 


Equations obtained from (4.15) for the indic5ated set 
of i and 3 are termed as adjoint equations . 

Making use of equations (4.8), (4.15) to (4.- 19), 6e 
can be written as, 


1 

I 

1=1 




1 


'1 k=1 


^ 8x!cW “ ^15^0)) 6X3^5 (0) (4.21) 


1=1 


"15 
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And 

since 6g^ and 

bgp are assumed to be 

zero , 

5 

k=1 

O' 

^ s 

= 0 

(4.22) 

f 

^gp 



1 

k=i 


= 0 

(4.23) 

PROCEDURE; 




The 

Optimum for the 

objective function' e 

can now be 


approached numerically by using the method of steepest 
descent. In this method one starts with an arbitrary deci- 
sion Variable vector and approach the optimum by moving in 
the gradient directions, which are in the direction of 
minimum. Thus 'find 6 e such that 

®new ^ ^old 

or 6 e . < 0 (4.24) 

This condition (Ecjuation 4.24) which is also necessary 
condition for minimum, can be met for (4 , 21),, if 

^ ( 4 . 25 ) 

6X^^(0) = - w V ■5x^|(o) - (4^26) 

where w is a positive weight function. 
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Then the potential computational scheme will he as 
follows: 

1 » lor assumed values of various variables which satisfy 
the system constraints, calculate coefficients of 6W^ and 
6Xij:j(0) for eqiiation (4.21), by solving equations (4*4) and 
(4.15) together with conditions (4.8), (4.9), (4.18) to (4.20), 
(4.22) and (4.23). 

2. If constraints are violated, decrease the value of w 
such that w >0 and all constraints are satisfied, 

3. How define ^^^^(O) such that, 

"llnew = ''llold - ^ 

^15^‘^^inew ^1 5 ^ ^ 1 0 Id ■ ^ ^ ( 0) ” ^15^°^^ 

Repeat steps 1,2 and 3 for new values of and x^^Co) 
imtill no further improvement is possible in e , 

It may appear that the coefficients of 6W^ and 6X^^(0) 
in equations (4.2l) can he calculated by integrating the 
equations (4*4) e-nd (4.15), end making use of equations (4.22) 
and (4,23) together with initial condition and boundary condi- 
tion (4.18) through (4.20). But, ri.j and ng in, equation (4. 19) 
cannot be calculated directly by any of the above equations. 

As in general cases, one can decide to assume certain values 
for n.j and TI2 check whether relations (4*22) and (4.23) 
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are satisfied after integrating (4.15) and (4.4). In this 
case one may have to search over the whole range of and 
^2 satisfy (4.22) and (4.23) , "because the system of 
ecjuations do not give any scheme to improve upon the assumed 
values. The problem is simplified hy using Green's function 
and Green's identity (15) for both Adjoint e q.uat ions and , 
variational system equations. 


Defining Green's function 
(4.^5') by, 


^ijk Adjoint equations 




df 


Ai 


ijl 


i 1,2f»» 
j and k .= 1 ,2, . ,5 


(4.27) 

with initial conditions 

^ijk^^^ ~ ^jk 

i.e. " ’’ ^ 

= 0 j / k 

The Green's identity for set of eqiations (4*27) is, 


x„(0) 


k=1 


^ijk ^^i^ 


^ik^^i) 


(4^28) 


Prom (4*19) one has. 
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i 

♦ 

and from (4»28), 


8gi 

8g2 

ds 


«N1 


Sg-, 

m 

3e 

• 

4r 

« 

8^2 I 
• 

• 


Sg2 

8e 

9%5 

»%5 

®%5 





^N11 

^N12 ••* ^N 15 


J 

T- 



^N2^^^ 

m 

= 

^N21 

, . ■ - . ^is [25 

% I 

- 

• 

• 

1^5(0) 

i 

j 

« 

• 

^N 51 

• 

* 

.... ^N 55 


* 

' 




Now making use of (4.30) and (4.29), (4. 20) gives, 


(4^30) 


^N-1,1^%-1^ 


%-] (0) 

^N-1 ,2^%-1^ 

» 

= 

♦ 

* 


« 

Ah4°>] 


^11 

■*■12 

= 

A21 

• 

• 


• 

^51 

• 

^52 
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and it is readily seen that can be expressed in terms 
of coefficients of, first matrix of right hand side of (4.29), 
first matrix of right hand side of (4.30) and second 'and 
third matrices of right hand side of (4.31). 

With the above procedure other values of Aj,^(0) and 
Xi;^(W^) can be known in terms of and ^2 and thus coeffi- 
cients of 5Xj^^(0) and in (4.19) can be known in terms of 
ri^ and rig. 


Now to calculate values of 6Xjjj(¥jj) in terms of 

and Tip which will be used to evaluate and rio from 
^ , ■ 1 2 

equations (4*22) and (4.23), one writes Green’s functions 
for the Variational System Equation; 


’"iok 



df 


il 


idi ax,. 


with boundary conditions 


r 


ijk 


(¥.) = 1 
■= 0 


D = k 


i = 1,2,...N ’(4.32) 

3 and k = 1,2, ,.5 


Greens identity for (4.32) is, 


6X..(W.) = J r. .^(0) 6X.3^(0) 


ID' 1 


^ik^ 


(4.35) 


And since there is also a 6W^ change in weight of catalyst 


¥ 


1’ 
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"I 

i 

! 


f r 

4112 **• Mi5 

■i 6X^^(0)' 

! 


f 

i . 

^121 

• 

• • 

• • 

) 

1 6X^2(0) 

1 i *• i 

' 



i 

j 

• 

• 

« 

: -i 

j 

f 




^141 

’’142 

^145 

' |5X^5(o) 

r 

! 

' t 

f„.6W^ 

J 



®12 

^13 ’ 

! 

i ^ 

®21 

• 

♦ ^ 

I 

1 

i •'2 

1 

i 

®41 

®42 

®43 i 

i 

. ! 1 


Here again coefficients B. . can be easily calculated.- 

1 J 

Making use of (4.9) and follovd.ng tbe above procedure 
6Xj^^(Wf^) can be calculated in terms of and ri 2 . And 

now substituting the values of in (4.22) and ( 4 . 25 ), 

. and can be calculated. 
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RESUXTS ilFD DISCUSSION 


In Tatle I, ■fcho outiputs froiti "bh© raodel ar© coiapaned 
with those from the two industrial plants. The details of 
each plant e.g., weight of catalyst in each reactor, feed 
and recycle gas composition, feed rate, recycle ratio, and 
temperatures and pressure of inlet streams are given in 
Appendices A and R. Though not all the desired data is 
available for the industrial reactors, comparison was made 
with two types of industrial units. One unit has three 
reactors in series and the other has four reactors in series. 
The amount of aromatics and paraffins in stream leaving the 
reactor were predicted with in 7.7' percent and '4.7' percent 
respectively for the two cases. The temperature of exit 
streams were also predicted by the model within , 4,8 percent. 
Production rate of hydrogen and hydrocracked products were 
predicted with in 2.1 percent. The prediction of naphthenes 
conversi'.'n was not satisfactory. A vain attempt was made 
to increase the accuracy of the results, by changing the rate 
constants appearing in equations (?,. 5) through (3.8) « However 
as the significant characteristics of the reactor systems 
have been predicted satisfactorily, it can be assumed. that 
the proposed model represents very well the true behaviour 
of the plants. 
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IPlBIiB I: 


SIMUIiAIIOlT RBSUIiT. COMPARISON OP 
AND PREDICTED VALDES. 


INDCSTRIAL 


and Ij 
P 
E 


*= Industrial 
*5 Predicted 
- Error ^ 


factor StreM tooma- Uaph- Para- H, pro- 
Se-t-np leaving ties thenea If ins duoed 

reaetor lb m/hr lb m/hr lb m/hr lb m/hr 


Hydro- 

cra#ked Temperature 
product Otj 



1 






xb m/hr 

Air 

1255. 


F 






1288*9 



E 






1.3 

Three 


^A 




. 


1320. 

in 

Series 

2 

P 

B 


- 




1330. 

.8 




170. 

8.9 

76.5 



1400. 


5 

P 

157. 

2.5 

77.4 



1374.5 



E 

7.7 

73.0 

1.2 



f.8 




74.7 

80.0 

76.0 



1271. 


1 

P 

73.9 

79.31 

77. 



1281. 

- 


E 

1.1 

.9 

1.3 



.8 



Ij 

109. 

39.6 

79.2 



1318, 


2 

P 

108. 

38. 

82. 



1322, 

Pour 


E 

.7 

4,0 

3.5 



.3 

in 

Series 


% 

129.4 

12.7 

75.0 



1368. 


3 

P 

126. 

11.55 

78.5 



1380. 



E 

2.6 

10.0 

4.7 



.9 



Ib 

140.1 

4.9 

42.5 

252. 

122.5 

1409. 


' 4 

P 

137.6 

5.07 

45.6 

244.6 

119.9 

1406. 



B 

1.7 

21.8 

.73 

2.1 

2.1 

.2 
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The profiles of conversions to aromatics, hydrogen 
and hydrocracked products and temperature for the plant of 
Appendix A are shown in Figures 2,3 and 4. It can he seen 
that in reactor 1 and 2, dehydrogenation of naphthenes to 
aromatics is dominant (Figures 2 and 3) and since this 
reaction is endothermic, temperature along the reactor 
length goes down. It can also he seen from Figure 3 that 
in reactor -2 more of hydrocracked product is formed than 
in reactor -1 hecause of high operating temperature in the 
reactor-2. Production of aromatics in the reactor -2 is 
less than in the reactor -1 , hecause dehydrogenation of 
naphthenes to aromatics is reaching equilibrium, even though 
the temperature of inlet stream and average operating tem- 
pemture are higher for reactor -2. 

From Figure 4, it can he seen that the conversion of 
hydrocracked product is dominant and that production of 
aromatics in very slow in reactor -3. One might expect 
the temperature to go up and concentration of hydrogen to 
go dOTO. But the results are contrary. The reaction (2.2), 
which is endothermic and consumes hydrogen in the forward 
direction, actually takes place in the backward" direction 
in this reactor. For example extent of reaction (2.2) at 
the entrance of reactor -3 is .000526 and at the exit it 
is -.008788. 
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^ Sh© siinulation ihocIgI was used "to s'tudy i^be effect of 
recycle rafio, and pressure ou fhe production of aromatics* 
The results are presented in Figures 5 and 6. Illie effect 
of the inlet temperatures of the streams on the production 
of aromatics j hydrogen and hydrocracked product were also 
studied. The results are given in Table II and Figures 7*8 
and 9* Additional plant data is taken from Appendix A. 

(i) Effect of Recycle Batio ; 

It can be seen from Figure 5 that as the recycle ratio 
rises, from 3 to 8, the production of aromatics also rises 
from ,3336 to .4018 mole of aromatics per mole of naphtha, 
and as it increases from 8 to 10 tbe production of aromatics 
remains almost steady. Since neaetion (2, l) produces hydrogen 
in the forward direction, one can expect that the aromatic 
production rate will go down with the rise in recycle ratio, 
as the increase in its value increases tbe partial pressure 
of hydrogen. But the results are contrary. As the hydrogen 
partial pressure in the reaction mixture rises because of 
rise in the recycle ratio, the exothermic reaction (2*2) 
bo comes fast. The temperature rise because of more heat 
produced in reaction (2.2) is sufficient enough to off set 
the decrease in the reaction rate of reaction (2.1). It can 
be further noted from Figure 5 that as recycle ratio rises, 
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• 2 mm II S EP^ECI OE laaPERATlIRE ON CONVERSION OE 
AROMATICS, HYDROGEN INN HYOROCRACEED 
PRODUCTS, 


T 

1 

Convesc- 
sion in 
Reactor 

1 

®2 

Conver- 
sion in 
Reactor 
2 


Conversion in 

Reactor 

“ 3 for 


"3= 

1285. 

T = 

3 

1320. 

T - 
3 

1355_. 

T^= 

i 3m._ ^ 
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' * 
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m Conver- m Conver- Conversion in Reactor 3 for 

1 Sion in 2 sion in -ifr-z nri nr^:: — nr zT ' ^ ' 
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.3182 

.3469 

.3911 

.4598 


1320. 

.5724 

,7216 

.8117 

.8748 

. 8968 

,8774 

.8479 



,0284. 

.0466 

.0675 

.1155 

.2354 

.5356 

1.0290 



.2334 

.2723 

.3038 

.3317 

.3637 

.4198 

.4818 

.1510 
1355. .3765 

1355. 

,6638 

.7656 

.8438 

.8908 

.8981 

,8798 

.8519 


.0429 

.0649 

.0925 

.1597 

.3333 

.7124 

1.1440 

.0175 


.2686 

.2900 

^3183 

.3454 

.3857 

.4503 

.4943 


1390. 

.7485 

.7953 

.8588 

.8888 

.8874 

.8772 

.8561 



.0732 

.1008 

.1391 

.2386 

.4865 

.9077 

1.2143 



.3024 

.3058 

.3313 

.3606 

.4105 

.4726 

.4961 


1425. 

.8096 

.7991 

.8477 

.8648 

.8632 

.8621 

.8430 



.1369 

.1713 

.2246 

.3666 

.6779 

1.0706 

1.2461 



.3343 

.3231 

.3468 

. 3779 

.4298 

,4802 

.4907 


1460. 

.8238 

.7686 

.8072 

.8189 

,8212 

.8246 

.8068 



. 2782 

.3162 

.3796 

.5461 

.8670 

1,1766 

1.2685 



.2022 

.2541 

. 2880 

.3178 

.5467 

,3927 

.4611 


1285. 

.5700 

.7126 

.8006 

.8610 

.8800 

.8600 

.8228 



.0405 

.0592 

.0810 

.1316 

.2591 

.5734 

1.0607 



.2291 

.2691 

.5007 

.3286 

.3604 

.4159 

.4785 


1320. 

.6455 

.7507 

.8295 

.8772 

.8851 

,8664 

.8380 



.0479 

.0697 

.0968 

.1627 

.5300 

.7087 

1.1454 



.2596 

.2849 

»3137 

. 3408 

.3792 

.4429 

.4914 


1355. 

.7253 

.7843 

.8509 

.8842 

.8838 

.8715 

.8496 

.1738 


.0647 

.0911 

.1269 

.2189 

.4512 

.8704 

1.2068 

1390. .4883 

.0361 


,2909 

. 3000 

• 3263 

.3547 

.40-25. 

,4669 

.4968 

1390. 

.7960 

.8035 

.8566 

.8770 

.8746 

.8709- 

.8515 


.1015 

. 1 340 

.1826 

.3118 

.6103 

1.0214 

1.2372 



,3210 

.3152 

.3397 

.3709 

. 4242 

.4809 

.4961 


1425. 

,8394 

.7996 

.8414 

.8544 

.8550 

.8574 

.8387 


.1826 

.2208 

.2827 

.4467 

.7787 

1.1299 

1.2521 



.3534 

,3379 

.3608 

.3908 

.4396 

,4839 

.4913 


1460. 

.8377 

.7703 

.8075 

.8192 

,8228 

,8265 

.8096 


,3685 

.4051 

.4667 

.6271 

,9265 

1.1968 

1,2677 


Contd. . ,45 
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TSonver- 
Sion in 
Reactor 
- 1 


1425* 


.2170, 

.5890 

.0717 


1460. 


.2564 

.6611 

.1415 


tT Oonver- Conversion in Reacto’r 3 for 

2 sion in ' m' — ’ 'iTi — ' " 'P — ffl _ ■ ' 'Vir _ 

Reactor 3~ 3" 3“ 3 ~ ■■ 

_...2 1285o 1320. 1^55. 1390. 1425. 1460. 


1285. 

.2344 

.6380 

.0763 

.2701 

.7296 

.0993 

,3006 

,8041 

.1288 

1320. 

.2572 
.7006 
. 0847 

.2817 

.7572 

.1113 

.3102 

,8225 

,1478 

1355. 

.2829 
c7653 
. 1049 

.2942 

.7797 

.1367 

.3206 
, 8340 
.1837 

1390. 

.3095 
.8180 
. 1520 

.3072 

o7894 

.1894 

,3321 

,8332 

,2492 

1425. 

.3370 

.8424 

.2618 

>3243 

,7821 

,3012 

, 3 1-81 
,8197 
.3672 

1460. 

.3751 

.8335 

.5096 

.3571 

.7611 

.5413 

"Mm 

, 5945 

1285. 

• c 0 

.6748 
c 1464 

.2819 

.7108 

.1749 

.5089 

.7695 

,2160 

1320. 

.2814 

.7214 

.1559 

.2901 

.7278 

.1383 

.3158 

.7790 

.2375 

1355. 

.‘3013 

.7669 

.1810 

.2996 
. 7402 
.2179 

.3242 

.7835 

.2774 

1390. 

.3229 

.7990 

.2440 

.3121 

,7438 

.2839 

.3358 
, 781 3 
.3510 

1425. 

.5510 

.8092 

.3959 

.3348 

.7397 

.4323 

.3574 

.7760 

.4943 

1460. 

.3991 

.0061 

.7039 

.3791 

.7315 

,7283 

,3992 

,7704 

o7690 


,3278 

,3612 

.4198 

.4779 

. 8464 

,8505 

.8335 

.8079 

.2021 

.3916 

.7848 

1.1883 

,3371 

. 3760 

, 4400 

.4868 

<■8543 

,8533 

,8417 

.8201 

,2419 

,4793 

.8900 

1.2247 

,3488 

. 3954 

.4605 

.4925 

.8553 

.8527 

.8481 

.8287 

,3087 

,6014 

1.0171 

1.2459 

.3629 

.4160 

.4758 

.4942 

.8471 

, 8467 

.8481 

.8292 

,4082 

,7404 

1.1127 

1.2554 

,3799 

,4320 

.4831 

.4930 

,3309 

,8337 

.8376 

.8198 

.5400 

.8674 

1.1746 

1,2625 

c:4056 

, * ,'’ ,79 

.4850 

, 4902 

.8117 

.8157 

,8186 

.8032 

.7328 

.9898 

1.2161 

1.2711 

,3358 

.3181 

.4423 

.4799 

.7947 

,7921 

.7844 

.7643 

,3244 

.5887 

1.0037 

1.2686 

,3441 

,3919 

.4554 

.4833 

,7978 

,7953 

.7923 

,7731 

.3689 

,6698 

1.0733 

1.2759 

. 3548 

.4077 

.4672 

.4854 

,7969 

,7964 

.7979 

.7792 

. 436 1 

,7673 

1.1379 

1.2791 

.3682 

,4225 

.4750 

.4860 

,7921 

.7944 

.7984 

.7803 

.5272 

,8639 

1.1854 

1 . 2809 

.3873 

.4359 

,4791 

.4859 

.7873 

.7911 

.7951 

.7786 

.6554 

,9527 

1.2157 

1.2825 

.4214 

.4544 

.4825 

.4856 

,7849 

,7891 

.7905 

.7770 

.8746 

1.0723 

1.2443 

1.2837 


Contd. . .46 
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How to read the tabl^ i 

' Upper, middle and lower numerical values of the three 
figures appearing together in the columns of 'conversion in 
reactor' stand for aromatics, hydrogen and hydrocracked 
product, produced per mole of naphtha respectively.' For 
example, aromatics, hydrogen and hydrocracked product pro- 
duced per mole of naphtha in the stream leaving, reactor 
one for = 1285 are .0559, .1644 and .0032, reactor two 
for = 1285, ^2 = *1437, .4191 and .Oi29, and 

reactor three for = 1285, ^2 “ "^520 and = 1355 are 
.2937, .8282 and .0718. 








■■■i 









W$M 




mm 

IMI§! 


' ■C '»**“ 


iHi 


{'f f r') lO ^'4 ‘.D CnI 
i-i- ii-; t 


(1i.:?U|,,?3 S&iiViA.!/"£ IO|,>u»j 
;( pai,np(,} .II.J u?s5o]|;AM I '-Asaiow 





50 


the increase in aromatics production rate decreases. At 
higher recycle ratio, the effect of extra heat produced 
because of reaction (2.2) becoming fast to make the reaction 
(2.1) fast, decreases. 

(ii) Effect of Pressure : 

It cap he seen from Eigure 6 that as pressure of inlet 
stream rises from 17 to 27 atmospheres, the production of 
aromatics rises from .5985 to .4131 and then decreases from 
.4131 to .3431 moles of aromatics/mole of naphtha as pres- 
sure rises from 27 to 47 atmospheres. Erom the Le Chatellier 's, 
principle one will expect the production of aromatics to go 
down with an increase in pressure, because both in naphthenes 
dehydrogenation to aromatics and paraffin dehydrogenation 
to naphthenes a rise in number of moles is resulted. As the 
pressure rises exothermic reaction (2,2) becomes fast and 
leads to a rise in the temperature of reaction mixture, and 
at high temperatures, the reaction (2.1) which is endother- 
mic becomes fast enou/.5:h to off set the decrease in the rate 
of dehydrogenation of naphthenes to aromatics because of 
increase in pressure. But as pressure rises, the rise in 
temperature because of (2.2) becoming fast is not sufficient 
enough to off-set the effect of the rise in pressure and as 
a result the production of aromatics goes down. 
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(iii) Effect of Temperature ; 

fhe effect of temperature on production of aromatics, 
hydrogen and hydrocracked product is given in fable II and 
figures 7> 8 and. 9. Figures 7, 8 and 9 represent a typical 
set of curves for one inlet set of temperature for reactor-1. 

It can be seen from Table II that the amotint of aromatics 
and hydrogen in the stream leaving reactor-1 is approxima- 
tely a linear function of temperature. The amount of 
aromatics produced in reactor-2 is also approximately a 
linear function of Tg for every v-^lue of T^. The amoimt of 
aromatics in the stream leaving the reactor-3 rises with 
rise in temperature of inlet stream to any reactor except 
for T^ > 1390, T 2 > 1355 and 1425°R. 

For the above excepted value a rise in inlet tempera- 
tures of all reactors, produces only marginal decrease in 
production of aromatics. 

For T^ = 1425 and T 2 > 1390, and T^ = 1460 and T 2 > 1355 
amount of aromatics is less in the stream leaving reactor-3 
than the stream leaving reactor-2 for T^ = 1285. At the above 
values of T^ and Tg, conversion of naphthenes to aromatics is 
very fast in reactor one a.nd two. The concentration of naphthenes 
in reaction mixture of renctor~3 for T^ = 1285 is low enough to 
reverse the direction abaction (2.1) as this reaction is 
endothermic. Since reaction (2.2) is exothermic, at low tempera- 
ture hydrogenation of naphthenes to paraffins is favoured 
and as n nf naphthenes is lowered. 

IJ.T. KANPUR 
I CENTRAL LIBRARY 

I 
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It is seen fron Fig. 7 that aiaount of aromatics 
produced at = 1460 and = 1460 is approximately inde- 
pendent of Tg* 

AS hydrogen takes part in all the reactions, its 
production does not follow any pattern (See Figure 8). 

The production of hydrocracked products rises with 
the rise in temperature of the inlet streams and it is seen 
from Figure 9 that it's amount is approximately independent 
of for T^ = 1460 and T^ between 1425 and 1460. 

OPTIMIZATION : 

The procedure outlined in Chapter 17 to minimize an 
objective function, e , with constraints on temperature and 
concentrations of the stream leaving the last reactor of 
the series was used to minimize the total weight of catalyst 
and operating cost for three, four and five reactors in 
series. The respective objective function for total weight 
of catalyst and operating cost, and the constraints are given 
in Appendix D. The results obtained by computer programme 
(sec Appendix E) developed in the present study, are given 
in Table III and 17 (For additional data, please see 
Appendix il) # 

When minimizing the total weight of catalyst and the 
operating cost for five ren,ctors in series catalyst weight 
for one of these reactors come out to be zero. Though tbe 
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TABLE , III : RESULTS FOR MINDIIZATIOS' OF WEIGHT 
OF catalyst 

SS = Starting solution 
FS = Final solution 


Reactor Total wt. 

3 4 -5 of cata- 

Ivst lb. 



SS 

?/. 

1 

1500.0 

1940.0 

Throe 

Reactors 


T. 

1 

1378.8 

1389.6 


FS 

w. 

327.9 

658.7 



L 

1402.3 

1380.6 


SS 

w. 

1500.0 

1940.0 

Four 



1378.8 

1389.6 

Reactors 

FS 

w. 

1 

1215.9 

1283,4 



'I. 

1 

1402.9 

1433,1 


ss' 

h 

1500.0 

1940.0 

Rive 


h 

1378.8 

1389.6 

Reactors 

FS 

»i 

0.0 

472,8 




1366.0 

1366.0 


13900.0 

1398.6 



17340.0 

12437.1 

1439.5 



13423.7 

6900.0 

1398.6 

6900.0 

1375.0 


17340.0 

831.4 

1460.0 

1387.0 

1408.7 


4717.7 

4633.3 

1398.6 

4633.3 4633.3 

1375.0 1374.0 

17340.0 

3241.7 

1391.5 

3242.7 

1438.9 

3263.8 

1372.6 

10221.0 
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TABLE lY; RESULTS EOR MINIMIZATIOH OE OPERATUffG GOST 

SS = Sorting solution 
?S = Pinal solution 







Reactor 



Operating 
cost in 




1 

2 

3 

4 

.5 

lakhs of 

Rs. per year 

Throe 

SS 

w. 

1 

1500.0 

1940.0 

13900.0 



32.36 


T. 

1 

1378.8 

1389.6 

1398.6 



Reactor 




w. 

1159.7 

1459.8 

12858.2 





PS 

1 

1367.7 





29.44 



L 

1391.7 

1410.3 






''■’i 

1500,0 

1940.0 

6900.0 

6900.0 




SS 






32.22 

Pour 


T. 

1 

1378.8 

1389.6 

1398.6 

1375.0 



Reactor 


W. 

1967.-8 

2036.2 

1445-3 

2336.6 




PS 

X 






17.75 



L 

1355.2 

1373.0 

1406.7 

1429.3 





'■■'i 

1500.0 

1940.0 

4633- 3 

4633 -.3 

4633.3 



SS 






32.39 

Pive 


Ti 

1378.8 

1389.6 

. 1398.6 

1375.0 

1374.0 


Reactor 


W. 

2298.7 

2462 . 3 

2302.1 

0.0 

1480.4 



PS 

1 






19.03 



u 

1357.0 

1375.0 

1409-9 

1360.6 

1436.7 
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■bh 2 ?e 6 startling solutions which wore tried led to one of 

going to zero, there can he an initial solution which 
can avoid such a situation. 

Since the iterations were stopped when one of the 
values went to zero, the final solutions given in Tables 
III and lY for the five reactors in series do not represent 
the real optinuin. When these values were used as a starting 
solution for the four reactors in series, approximately the 
same results were obtained as given in Tables III and lY for 
four reactor cases. 

From the results obtained it is concluded that one 
should use four reactors in series. But the result can be 
different if one maximizes for venture profit. In operating 
cost the number of inter heaters, which depend on the number 
of reactors being used is not considered but will figure in 
venture profit. 

During computation, it was found necessary to reduce 
the scale of temperature and naphtha feed rate by factors of 
10-40 and 50-100 respectively. It was also found necessary "feo 
reduce the scale of optimization function by a factor of 100 
when minimizing the operating cost. These scaling were made 
necessary for, in the development of algorithm, second and 
higher order terms have been assumed to be negligible as 
compared to the first order terms. In order to reduce the total 
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nuinlD6i' of i'toica.'tioii s if. i/vas found usoful "bo us© a. lowei? 
factor for temperature and a higher for naphtha feed rate 
initially and then keep changing these factors after some 
iterations. 

The values of constraints g^ and- g 2 were .046 and .5731 
respectively at the starting of the iterations for optimization 
in all the cases considered for the optimization study. These 
eq.uallty constraints were allowed only a change of + .0005 and 
+ .005 respectively from the starting value. If the change in 
the values of g^ and gg was more than the above given value, 
w*, the multiplication factor to determine the magnitude of 
changes in decision variables, which jhad an initial value 4, 
was reduced to half each time till the constraints met the 
prescribed limits. 

The iterations were stopped vdien changes in objective 
functions for the weight of catalyst and the operating cost 
were less than 1.0 and 50.0 respectively. Por both the 
objective functions, the time taken on IBM 7044 to reach the 
final solution from the starting solution were approximately 
20, 35 and 15 minutes respectively. 



GELAPTER TI 


CONGLUSIOHS MD REGQMMEITDATIOITS 

Tho 1118.13116111813108,1 luodel developed, foi* Ihe catalytic— 
reforming process predicted very well the temperature and 
conoentratioi-i of various components in the exit stream of 
the reactors, for two different industrial units. 

Simulated results show that as recycle ratio rises, 
the production of aromatics also increases. Maximum produc- 
tion of aromatics occur at a pressure of 27 atmospheres. 

The rise in temperature of inlet stream increases the pro- 
duction of aromatics and hydrocracked products. Erom the 
simulated results, it is recommended that for the plant of 
Appendix A, the temperatures of inlet streams should he 
changed to T^ = 1390, Tg = ^590 and = 1460 from 

the existing values, to increase the production rate of 
aromatics to 135.2 moles/hr from its present value of 112.2. 

The increase in temperature will lead to more load on inter 
heaters* So one has to make an economic evaluation of the 
process before taking the final decision. 

If one wishes to use the minimum weight of catalyst, 
or to have the minimum operating cost to produce aromatics 
at a rate, being produced in plant of Appendix A from 
naphtha (its composition and feed rate are given in Appendix A) 
one should use four reactors in series. 
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In the existing plant (Appendix A) one can use less 
weight of catalyst or operate at lower operating cost for 
the Same aromatic production rate and the same q_uality of 
reformate. The weight of catalyst can he reduced hy 19.8 per 
cent (See Table III) and operating cost hy 9.05 per cent hy 
changing, the catalyst distribution in the reactors and tem- 
perature of inlet streams to the- reactors. 

KSCQMmTtilTIOhS FOR FURTHER STUDIES : 

1» -^n the present study, inter heaters have not been 

included in the simulation. A mathematical model representing 
the behaviour of inter heaters can he combined with the pre- 
sent simulation model. 

2. Mathematical model representing the flash unit can 
also be incorporated with the present simulation model. 

3. Re cycle ratio and operating pressure can also be 
used ao variables in the optimisation function. 

4* for design optimization of Platformer unit, a venture 
profit concept can be used as an objective function. 
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*;PPENDII A 


PIulIT DiJJ A FOR gHREE RliklGTORS SET-UP 

Food = 272.2 raolos/hr, 

Rocyolc; ratio = 7.95 


Reactor 

uuacia^x Qfiijalyet ib 

Temp, hi inXe't 
stream OR 

treasure' dS 
inlet stream 
atm. 

1 1500 

1378,8 


22.0 

2 1940 

1389.6 


20.5 

3 13900 

1398.6 


18.6 

Naphtha Feed Coniposition; 



CY“Ptirtaffina 


2.2 ^ 

% by wt. 

CrpNaphthcnoo 


2.9 

ff 

0 ..^-Aromatics 


2.0 

ft 

Oj^j-Paraifins 


29.0 

ft 

Cg-Hc.phthonuo 


29.3 

f! 

Cj,^-Aroinr,itico 


12.6 

!! 

C(.-Faraffins 

V 


18.4 

fl 



3.6 

!f 
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Rocycle Gas Composition; 


Hydro gon 

86.0 

Mole 

Methane 

6.0 

n 

Hthane 

2.8 

If 

Propane 

2.1 

ff 

Butane 

1.0 

ff 

Pcjitano 

0.5 

tr 

C^^-Paraft‘ins 

0.6 

»f 

Cg~Cp Hydrocarbons 

1.0 

ff 



.IPPENDIX B 


PLi^T BOR POUR RK.CTORS SET-IP 

Feed = 232.8 r.ioles/hr 

Recycle ratio =6.9 


H(}ru^toi* 

Vfcight of 
('at.'ilyat lb 

Temperature of 
inlot stream °R 

Pressure of 
inlet stream 
atm. 

1 

1180 

1403 

39.9 

2 

1180 

1403 

39.5 

3 

. 3230 

1403 

39.1 

4 

4150 

1403 

38.1 

Naphtha 

feed composition; 



Cg-ilaraff ins 

0.3 °/ 

0 hy wt. 


Cg-Naphthenes 

0.7 

tf 


Cg-iiroinatics 

0. 1 

fl 


Cr^-'Paraffins 

5.9 

ft 


C^-Naphthenes 

9.3 

T! 


G^-Aromatics 

3.7 

tT 


CQ-Paraffins 

23.9 

tf 


Cg-Kaphthenes 

26.1 

ft 


Cg-ilroinatics 

11.5 

ft 


C ^-Paraffins 

15.0 

tf 


Cg-Naphthenes 

1.9 

ft 


Cg-ilromatics 

1.6 

tf 


Recycle G-as Composition; 

Hydrogen 86.4 Mole % 

Methane 5.1 ” 

Ethane 3.6 *' 

Propane 2.1 " 

Butane 1 , 406 ” 

Pentane 0.666 " 

Hexane and higher 0.698 " 
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HELII CAPAC PITER 

let It be ass«nea that SPi. + SP2..T, SH1 + SUa’'!’ 

SH1 . SH 2 .I are- specific heats and; XP., a’ 

X-,. XHI aie .Ole fractions of, pa.alflps, naphthehes^Po:atlos 
abd hydbosen in the feed to fl.st beactor whefe .f is pa^ber 

of carbon atoms. Xnd sin^ft ^ , 

and higher compounds are 

hegllglble, the specific heat for various classes of oo.pouuds 

^ .reactor e.g. , Cg and higher paraffins, 

P thenes, aromatics and and -lower paraffins, hydrogen and 
hydrocarbohs produced fro„ hydrocraClpg, car be wrltteh as: 


2 SP1.,XP. 
1 1 


9 


I 

XR. 

0=6 

. 1 

9 


2 

an. .xin 

i=6 

1 1 

9 


2 . 

m. 

i=6 

■ 1 

9 


2 

SA1 .Zil. 

i=6 

1 1 


2 SP2. .XP. 
i=6 ^ 1 


9- 


2 

XP. 

i=6 

X 

9 


2 

SN2. .XN. 

i=6 

1 1 

9 


'■■•2 " 

XN. 

i.=6 

1 

9 


2 

Sii2. . XA . 

i=6 

r\ 

1 1 


A Xil ' 
i=6 ^ 


A m. 
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^pHCI 


5 

s 

i=1 

5 

2 

i=1 


SPl,.iLP, 
1 1 


2P . 

X 


+ 


5 

2 SP2,*XP. 
— ■1 i i 


5 

2 XP. 
i=1 ^ 




= SH1 + SH2.T 

-I 5 ^ 

= i H 


5 

2 SP2..3) 

i=1 ^ 


Values of SP1^, 312^, SAI^^, SA2^, SH1 and SH2 for the 
catalytic reforning temperature conditions are given helow (l) 


SH1 = 4.02 

SH2 = .00024 


i 

SP1. 

SP2^ 

. 3171^ 

SII2i 


S42i 

1 

5.34 

.00678 





2 

8.97 

.0117 





3 

13.30 

.016 





4 

18.77 

.0222 





5. 

23.53 

.02492 





6. 

28.28 

.0292 

21.48 

.02958 

18.55 

.0184 

7 

32.98 

.0326 

25.62 

.0344 

23.57 

.023 

8 

37.86 

.038 

30-50 

.0388 

28.18 

.027 

9 

42.53 

.0424 

35.07 

.0432 

33.5 

.0311 


References: 
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American Petroleum Institute, Research Project 44 




APPENDIX D 


OBJECT I7E PUITOTIONS AND CONSTRAINTS 

1. Objective function for total weight of catalyst: 

= + ... % 

2. Objective function for operating cost: 

{Yft, of catalyst (lb)) (Cost of catalyst (Rs/lb)) 
life of catalyst (Yrs.) 

+ heat duty (Btu/yr)- heating cost (Rs./Btu) 

+ Eg,, pumping cost (Rs./(lr. lb mole /hr.)) 

Cost of catalyst (l) = 140 Rs/lb. 

Life of catalyst (i) = .9 years 

Pumping cost (2) 

= 38.2 ((pressure of stream entering the first 
reactor of the series/pressure of 
stream leaving the last reactor of 
the series) - 1.0) Rs/(yr. lb m/hr) 

Heating cost (2) = 390.0 Rs./(million Btu) 

3. Constraints: 

1. (Wj,) 

2. g^ = Concentration of aromatics in the reformate 

leaving flash unit 

Moles of aromatics in reformate 


Moles of reformate 
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It can be seen from the recycle gas composition that 
moles of and, higher hydrocarbons in recycle gas are appro 
ximately ,0134 moles per mole of hydrogen. Assuming that 
whole of the hydrogen is seperated in the flash unit, it can 
be said that .0134 mole of hydrocarbons/mole of hydrogen 
produced, are seperated in flash unit along with hydrogen. 
Thus, 

Moles of reformate per mole = XAI + MI + XPI 
“ “ ^14^^%^ -.0134 (3Xj^^(l^) - 

§2 = (XAI + Xj,,,(5rj,))/(XAI+XHI+XPI-.0134 . 3X^,|(W^) 

+ .0134 Xu2{l%) + (.0134 ^ - ^)X]J3(W^) 

+ (.0134 f - )) 

References : 

1. Nelson, W.N. , Petroleum Refinery Engineering, 
McGraw-Hill (1958). 

2. AlChE Student Contest Problem (l9i'^) . 
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4FFENDIX e 


jrj2 cgg:i;',ti‘' " t,^ag s > ..name ashok kunar sharma 

IBJG6 
IB FTC 

SJ''HJL^T!T\ 

THIG PF ;!.T- -.v'-. ; r j n pp dci ',h- OUT PUT STRFAMS FROM EACH REACTOR FOR 
ClViM A' CYCL“ C '.S FLCK PATT , hAPHTHA F~ED RATS, TEMPERATURES AND PRcSSURES 
FCP !''L ;:Tp ■ U’ TC T.Ch R[ ACTOR AND WEIGHT CF CATALYST IN EACH REACTOR 

JR IS ‘!U?'3 P OF R' ACTORS IM S PIES 

p, p A.Mp o-rv ^01”=^ CF NAPHTHA FE^ED AND RECYCL- GAS PER HOUR 

FP 15 PR.- SSL"' '’F TH"'- J iLLT STREAM TO VARICUS REACTORS 
TT ATD wC ■'^P'' T 'YF -RATURHS CF INLET STREAMS AND WEIGHT OF CATALYST IN 
■'ACH P .actor 

C6P1N- I;: f^CL., FHACTIOM rjF C6 AND HIGHER HYDROCARBONS IN RECYCLE GAS 
INC IS •rJ^'^PEP CF PARTS TJ WHICH A REACTOR IS DEVIDED FOR INTEGRATION 

DATA Cards 

MULT FhACTICNS UF R CYCLE GAS COMPONENTS STARTING FROM METHANE TO HEXANE 
AND HYCR'jGrN IN 1- FG.** FORMAT 

WEIGHT PTR CENT ACTS OF PARAFFINS (06 TC C9), NAPHTHENES ( C6 TO C9) AND 
AND ARnOATICSiC': TO 09) IN I F8«.4 FORMAT , 

NUMB'-T OF F^E ACTORS (I ft) 

TLMPFRATLR'' VSCTOP ,W’PIGHT OF CATALYST VECTCR AND PRESSURE VECTOR (SFl-.®' ) 
MOLE FkACTIC IS -'iF CS.AND HIGHER COMPONENTS IM RECYCLE GAS (i!‘F8o4) 

Ffcf.D AND R-. CY ( BFl' o ) 

REAL MWF,MUR‘-,r; 

DIMINSION A AA ( ■^) , UAH{7 ) ,AAHC (7) 

DIMONSirjN CF'\{:.),CPN(2) ,CPP (2 ) , CPH { 2 ) , CPHY (2 ) ,CPHYI (2) 

DIMENSION TT{4),FP{4) 

COMMON /3Li/ X1,X2,X;: ,X4,TI ,XAI,XPI,XNI,XHI,XHYI,AR 

COMMON /BLD/ IJ,TT,AB 

COMMON /BL V A i , A2 , A ^ , A^i , P , FT , N 

CCMMCN /BLR/ 40(2,=;) 

COMMON /BLft/ C{ - ) ,H 
COMMCN /RLC/ INC 
CIM)-KSir'N .sCC^) 

CALL CP”'V 
C { 1 ) =2 c 

C(2) =?.. 

c {?•)=:. 

FT=Fr:FL) + R'' CY 
r.^.AD e< , IP 

READ 6.,, (TT { I) ,WC{ I ),PP( I) , I=a,IR) 

ftRi FORMAT { IS ) 

£I FORMAT (BF].' .', ) . ' 

RHADAliFE'IStP.FCY 

PfiADft'.tCftFINR 

■:AI = YAF^-C0rJV+R:;:CY-»C6PINR/FT 


O AJ kJ 
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14 


IMC=6 

^"^=2 p'-. { I ) * r* 5 T + C'’? ( I S :p I +CP h { i )«/: MI+CPHC I y «XHI+CPHY I ( I ) »XHYI 
<\'>CPi ( I )• CF K I ) +'* •'JCPHC I ) 

«^P=CPP ( I )“'CF’!( I )*’CPH ( 1 ) 


!\ 3 =CPHY( ! ) ^N/:. 
A4=C PHY( I j 
AD(I )=*4. 

AD(I tP )=A. 

AD(I 

AD (I t^-0=A': 

A 0 ( I , * y=A^- 
CQNTINU,*; 

A 3 = FT/ 2 ‘'.'’c.? 

PP INTI 


■ cpps I )'• f ■■p- Ao)-’iCPH( I ) n> 

• CF ( I )-CPH{ I ) :^K/3o 


:-R' 


H 


FORMAT {////I 
1 PARAFFIN 
Xl=Oo 
X 2 =' « 

X3= c 
X4=l-o 

CO I 1=^.1 IP 
TI=TT{ I) 

H=WC (I )/FLCAT{ 
P=PP{I) 

CALL SIMULA 
AA = ( XA I +XI ) ^‘ FT 
AP = { XPI + X;.- X,J) ^FT 
AH = ( XHI + Fa'i^Xl- - ( 4- 

AN=( XNi- x;;- <. i )*-ft 
AHY= {XHYI + ( > 


?.CTCr-> v;T cat a 
HYCROC*//} 


T IN P. T 


OUT R AROMAT NAPHTH 


1C) 


)*X3/2c- N-»X4/3 o)*FT 
. , . j y-^FT 

PR INTI T» I,VvC(I)?TT{iy»TI»AA tAN ,APi AH, AHY 


13 FORMAT (1 
I CONTINUE. 
STOP 
FND 

IBFTC SUBl 

subrout Ii J£ 


X » I '’!• t F L 2 « j. j 


FCo ■ , F9ol » F7o i »FBo 1. » F9«lt F1»’8 1 1 F8o 1/ ) 


SIMULA 


this SUBROUTIU" I 


TO SIMULAT, THC RFACTORS 


REAL N,K14,5) 
DIMENSION PFK ") 
DIMENSION TT(4) 


COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 


/BLl / 
/3L2/ 
/3L3/ 
/BL4/ 
/BLA/ 


,X^+»TI ,XAI tXPI ,XNI,XHI, XHYI , A!i 


X 1 , X, 2 t >' 

J,TT,A3 

A 1 j A B f A ,1 , A 4 t P » FT ? N 


AD (2, 5) 

C { I ) , H 
/BLS/ INC 

;>'347:.. o/T) 
AZ 1( T ) =t:XF ( 4 *16 1 > •46' 4 /T ) 
A12( T) =‘.'XP (3Fo 93*’ED6'.i o /T) 
A22( T) =FXP ( SO 8 /T-'TsIZ ) 

A 33 t T) =LXP( 4 2 o « «/T) 


Ali( T) =lXP ( 2' 


7i 


/.NI = X''1F*F/ iC/FT 

\/p T - V p F c r /P T 

XHYI :XY*C6PUR/FT 

XHI=\H-',^^>?::CY/F"- 
:*l=(MwF-?c-.-'X )/lM^ 

F 1 ( p ) = -v p-4 ( '/I' p-^ ■> ’ •'»■ A'* X 2 / { ^:^ 2* A**3 ) ) / ( FT*A 5 

F2 ( P ) = ^ ’* ■'■P * ^ •^X'- .■■! h P/c^;- ) / ( FT*’ A ) 

F3 (A ) = ".■:.•> -X/ (F'^^'-a) 

F4( A )= ('^T-f'-A) 

F5 { C I ) = (• C;,-* P?. e+-.9 s -*02+ ( N‘’*Fo )*Sl^'Y o*C3’<-N*’22 

CO 1 I I = ; > IKC 

T=T1 

X1C=X' 

X2C=XP 
X3C=<3 
XAC=a 4 
DO 2 i=;u* 

A= ! o +3 * •**’X 1 C*" 1'IC‘ 

XA=XAI +XiC 
XP=XPI +'<PC-'X'X 

XH=XHI +3c ’C*-’ ( ‘A’ .>0 )’»’X3 C/Xo'’N/'3o*’X4C 

XN=X*a-- {X-C4XP.C + X<i!C) 

F1=A11 CT) 

'■:2 = A21 (T) 

F.3=A'i2(T5 
a4=A22. (T) 

E5=A31(T) 

Cl=Fl{P) 

C2=F2fP) 

C3=F3(A) 

C4=F4(A) 

Ai, =AD( ), t x)»T+AD{2f ) 

A i = A D { ', j 2 ) ■'* T ■FA 0 f 2 1 . ) 

A2=AD(,i, + t) 

A3=AQ(:u4 )«T+AD( 2 ,'f) 

A4=A D( 1 j f )’*‘’T+AD < 2 » P ) 

C P = A ' + A ] *• X C -F A 2 'J!- X 2 C -F A 3 *■ X 5 C + A 4 * X 4 C 

K(I»S)=C:i 

K(I,C)=C2 

K( It 3)=C3 

K( I, 4) =€•'! 

K( It 5) =F3(C1 ) 

IF( IeFQo4) GO TO •: 

XXC=XJ +H'*’K( 1 1'~ ) /C ( i } 

X2C=X2 + H^K( 1,2 )/C( I ) 

K3C=;r/+H»K( 1 1‘3 ) /C ( I) 

X4C=X*' + 1-I-»K{ I i4) /C ( I) 

T=TI ■fH’^K ( 1 1 f: ) /C { I ) 

CONTINUE 
DC 3 I ==1. t 5 

PHK I) = ( K (! 1 1) -Flo^KI 2t I 5 + 2® «K (3,1 )+K(4 tI ) ) /6» 

X1=X L+ O'* P H I ( 5 
X2=X2+h^PhI (2) 

X3=X3+ HtFpH I ( 3 ) 

X4=X 4+ H* PH I ( 4 ) 


-<,*C4/3<, )/CP 
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4 { /. , F-S 'F: .4, FI ol) 

cn'”!!!:] 

pr.:T!j?'l 

"■■'JO 

13FTC SUB? 

SUBRf'lUTI j,: C">"-V 

C0M‘l-r.) /BL.''';/ ,CPN{ 2) ,CFF{2) ,CPH{2),CPHY(2),CPHYI(2) 

CIMMIM /BL. ^/ ;:PF,XMF,XAF»XFYR,XH2Rt'^WF,MWRF 

0 1 M -■ -v S : n - I w C " { ‘ ) t W C P ( 4 ) , WC '-H ^ , W C A ( 4 ) , C P ( Z » 1, f : 3 1 C N ( 2 , 4 ) , C A ( 2 , 4 ) 

P;-:4L nr>(7 3 » '’CP( 4) ,-icm(4) ,PCA{4) 

AL M ’•< F > Y Vi P '■ 

OAT/. CP/, , -.Oe 72,2£c n,34B2-' T37^e,44o53»48o2^,54,33,59«B< , 

I fl ^'6, 7 P,. 5', ^40, Blot ), 9203. ,9So 93 »1* £«, 56 ,7o. 35, 7o« 7B/ 

DAT/' ,;;9o 2f^5.,.I2,79c3' »86*28,59« 57»97o35/ 

DATA C//-:.7',- 'o <1, ■ "9 , .-'So /2i62o2F»67oi5»72o7i .t78«3/ 

R . A D i' . ? U 1 . 

R £An 2 , WC F* , 'ft C ' ! , w C A 
2 FDRMATd' Fd \) 

A = ''.o 
R= o 

CO 8 i=:,,e 

MCR( I) =rtCP. ( 13/1 o 

A = A+*,4c 

a B=B+YCP { 1 3-»a 

MCP( 7) =WCR ( 7 3/ 

MWRF=P+MCR(7 5-»2c 
XHYR=1 ,-**YCR (T) 

XH2R=MCR{7) 

A A =7 2<j 
XAF = i , 

,>CNF= o 
XPF = ' t 
CO 1 1=.. ,4 
AA = A.^:,'-w, 

cc=aA“:4 

MCP( 13 =WCP( I )/ iA 
MCM{ I) =WCN{ I 3 /BB 
MCA{ I 3 =WCA { I ) /CC 
XPF = aPF+MCP { I) 
xnf=:';'JF4!'/cn c i 3 
1 XhF=XAF + wCA ( 1 3 

A=/YPF+XNF + X AF 
MWF=:0'";c./A 

xpf=::pf/> , — ^ 

XNF=;'CMF/A ■ 

XAF=XAF/A 
00 4 1=1,^ 

MCP( !3 =,‘1CP { I 3/A 

MCM( I3=YCK{I)/4 A 

4 MCA( I)=nCA(I5/A ' V 

4 = !''o 

B ' c 


6 3ol<' 7 


CD='-' o 
DO 5 

M=“( CPllil )-ZP i ■.,!)) /I ■ o 
C=CP( o 

D = D+M 
DD=DD+C 
ft = A+M«-;iCK ( I ) 

5 3=8+C--^MCP ( I ) 

CPHYI(1) = A/XHY- 
CPHYIC ' )=a/.XHYR 

CPCmI'. )‘'C?( I, , ))/:.■ 
CPH{ 1) =Y 

CPH( 2) =CP{ 1 , : )• M*' . 

CPHY{3)=C/?a 

CPHY(2}=DD/£. 

A = »',c 
R ” © 

C=i e 
AA=!':o 
8B='’'o 
C D ” : <S> 

DO 6 I =l 
J = I + 5 

H=«’(CP(1,J)~CP<2, J))/l r, 
C = CP (1, J)- , 

A = A+M»MCP(n 
AA=AA+C«MCP(I) 

M=-( CM (1,1 )-Ct! ('’,!)) /l^ o 
C = CN(1 ,I , 

B=B+M»HCM{ I ) 
eB=BB+C»^MCN ( I ) 

«=■“( CA (1, 1 )‘’CA (•:,!)) /I • c 
C = CA{1, 1 )-’N-97 ' o 
D=D+M»MCA{I) 

6 CD=nD+C»MC A ( I ) 
CPP(1)=A/XPF 
CPP( 2) =AA/XFF 
CPN( 1 ) =6/XNF 
CPN( 2) =BB/XNF 
CPA{1)=D/XAF 
CPA( 2) =DC/XSF 
CPP( 1) =CPP{1 

CPN( 1) =CPM { 1 )/Ic B 
CPA{ 1) =CPA{ 1 )/..ot5 
CPH( i) =CPH(1 )/..o 3 
CPHY(1)=CPFY(l)/;.c: 

CPHYK X )=C?hYI ( n /Ic c 

RETURN 

END 

EMTB.Y 
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JOB CGG: TIM T,P-Vo S j tNAME ASHCK KUMAR SHARMA 

IB JOB 

IB FTC MAIM 

PR'lGf-'A'M* r PC=> i^IA’IMlZATIOfi CF WEIGHT CF CATALYST 
IR I" 4J‘'‘P.;F! IF R^iCTHRS IN STRIES 

FT'D fi'iD P CY AFT 13 MDL'-S CF NAPHTHA FEED AND RECYCLE GAS PER HOUR 
PP IS PFI.S:lP^' of TH:: INLfcT STREAM TC VARIOUS REACTORS 
TT A'’,n vjC T'-.MPTR;, TURFS CF INLET STREAMS AND WEIGHT CF CATALYST IN 

EACH '^.-aCTCo 

INC IS NIP-B'£R IP r.'KTS IN WHICH A REACTOR IS DFVIDHD FOR INTEGRATION 
C6PINK IS '•'iCL': FRACTICN 'IF CH AND HIGHER HYDROCARBONS IN RECYCLE GAS 
ABl IS N'JYCFP IP MCLI 5 OF NAPHTHA WHICH ARE USED AS BASIS OF CALCULATIONS 
FOP OPT 2 ATI V-i 

ABC IS CCALI.IG FiCrOP FCft T EM PER ATURF 

VALUES UF AE . -iMC /..,3C AP”' CHANGFC CURING ITERATIONS AT SUITABLE PLACES 
THOUGH TH*-: GIV:M STA'^aiNG VALLES HAVE E5EN FOUND TO Be GCQC, BUT IT IS 
ADVIS"n TC TRY OTH R STARTING VALUFS FOR THHSb FACTORS TC MINIMIZE TOTAL 
MUMBTR OF IT R df.NP 

CIV IS A FACTOR Tr D-CID" ‘'AGivITUDF CF CHANGES IN DECISION VARIABLF.S 
SO THAT ALL Tl-r* CONSTPAINTS AOc SATISFIED 

ITR? IS MA'UwU/1 NLYS' R CF nFRATIONS FGR WHICH THE PROGRAMME IS ALLOWED 
TO EX'-CUT,: 

FR:. A'pn -R.-- - AMCU 'TS PF hRRCRS ALLOWED IN CONSTRAINTS CNE AND TWO ON 

THE FJi'lAL product ST'-FAM 

OR, IF TH CHA'IGC IN' OBJECT I VC FUNCTION IS LESS THAN ER PROGRAMME IS 
ThRMINAT-D AUTO m ICSLLY and THr SOLUTION VECTOR AT THIS PLACE IS PRINTeC 
AS OPTIMIZf-C RL CULTS 

IT IS ADVISED THAT MANY INITIAL SOLUTIONS EE TRIED TO ENSURE THE VALIDITY 
OF TH'^ P-SULTS 

DATA CARDS 

MOLP FHACTICs'S OF PCCYCLF GAS COMPONENTS STARTING FROM METHANE TO HEXANE 
AND HYCFNnGi,.N I j 1' F3-4 FORMAT 

WEIGHT RCRCE'lTAGfS CF PARAFFINS(C5 TC C9), NAPHTHENES ( C6 TO C9) AND 
AND AR iMATICSHCe. TC C9) IN liF8«4 FORMAT 
NUMBER OF REACT'IRS (lo) 

TEMPFRATUPR VCCTDR , W -I GHT OF CATALYST VECTOR AND PRESSURE VECTOR (BFlto* ) 
MOLE FRACTIONS OF CA AND HIGHER COMPONENTS IN RECYCLE GAS (ir'F8o4) 

FEFD AND r;;cy (3F:. c ) 

THFSC TT ttNC WC VALU'.'S ARE: LS:=0 AS A STARTING VECTOR FOR OPTIMIZATION 
at THrSr VALUl. S ALL IH- CONSTRAINTS SHOULD BF. SATISFIED 

CO'^MGN /BLl/ XlfX2fX0,X4 
COMMON /8L?/ A.<.tAj,A‘',A't,P,FT 
common /BL6/ C ( ■'0 ,H 
COMMON /BL7/ IJ,FL(S»5) 

COMMON /BL2/ I NC » I NC i 1 1 NC2 
COMMON /Blv/ ABC 
COMMON /RLCC/ N 

COMMON /3LU/ C0( A ) ,Cl,C2,Bl,e2,B3,DI V, IR, IRl 
COMMON /BL12/ T I , X A I , XP I , XN I ,XHI ,XHY I , AC 
COMMON /3l:.3/ T t n ( 5 ) , WI N ( 5 ) 

COMMON /9L\A/ WC,A 
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cnf^rm'i /3 l f/ ) iCpn(2) ,cpp(2),cph(2),cphy(2),cphyi{2) 

CGYf'iv'J /nL„c/ :PF, .:'JF,XAF, Xl-YR,XH2Rf MWF, NWRE 

cop?iDM /:a:?/ \r;{ 

CIX’:-\F r‘N ?F f i ) ,WC C l,TTI 55 

1 1 : , ^ , s \ a( ^ , ft (5 ) ,pw{5 ) 

CIW'v^FI^l 1 /;8 M ■ ) , 

PFt.L 

data PT,P'A/iH T. , '.il T2 , SH T3 ,6H T4 , 6H T5 ,6H Wi , 
16H wP ,.'F D j 'T! W4 , £H 5v5 / 

-F=j. 0 

CALL CP.Y 
"■Rl=o' 

■■R?. = c' ■ L 
WC 2- if. + ■ 

JC4 = 

■A3i= i,o :. 

Aec=? . 

CiV=A:, 

DIV1=DIV 
R:\hDs ,n 

p : AD A f ( T 7 ( T ) , rtC { 1 5 , PP ( I ) , 1 = ^ , I R ) 

6'' FORM«T(If) 

61 format { PF ^ ® 5 
IR..= IP-:: 

IT-=<, 

C(1)=Pv 
C( P) =2<r 
C(3}=!,c 

RcADf''l}F- 'DjKrCY 
R FAD 61 , CAP I NO 
AFT=FF:~D+R 7CY 
A=F:=-n 

CO 1 1= , IR 

18 wc(i ) = wcn )«A3 :/a 

PRI’^Tl 

3 FORMAT C H’j } 

Z FOP.MAT( ////I :,T2( 5 /.A Xi 3H^ ,7i X , IH*/1C*X , IH* t3uX , #FEF.D CATA*,3f 

1.x, XilF ’-,-' C , X,72( ,iH*) / 10 X,IH«, 7 FX, 1 H«/Ii X,1H*,7 
PRIMT1 7, F. ::C ,R. CY,MWF,MkNRr 

17 FORMATCj :<,1H*,'"'X,4,nAPHTHA Fi'£D LB HOLE/HR =« , F8o .1 , 17X , 1 H«/1LX , 1 H 
1^,70;, ;',»r=CYCL7 gas LS mole/hr =-»^,F8.1,l7X, 

2IH*,7' X, XtlH^tlTXt^LB MCLSCLLAR WT OF MAPHTHA =*, F£a 1, 17X, IH 

3*/;' K, lH*,7n <, X,1H*,12X,*LB MOLECULAR WT OF RECYCLL GAS =-»,F 

t, 17;:, .7, ) 

CD(1 ) = ;:AF'»FO-'D+r iCY4*C3PI'']R 

CD(? )=XXF#FErO 

DD(3 )=XPF#F1"D 

CD(4 )=::HYR'sR'''C Y‘>RtCY>C6PINR 

A=XU3R*!’Rf';CY 

PRPIT] 4,CC,,ft,PP{l) 

19 FOPMAT(’;,<' X,IH»,?3X,"»AR0MATICS LB MOLE/HR =#, F8« 1 , i7X , iH»/ir'X , iH* , 7 

X,!H-^?2X,*-NAPHTHeN£S LB MOLE/HR = «, F8o 1 , 17X1 H»/ir XlH^^7(„Xi, 
K:h-’!P7X*^p\RAFFI'7S lb mcLl/HR =*F8ol,I7X,lH*/lFX,lH«,7« X,lh*/1 
•Ji X,*HYi)ROCARBQNS LB MCLE/HR =■», FB® 1 , 17X, 1 H*/ 1* X , IH* ,7' X , 1 H 
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4#, /I X, :.h*, LB MCLS/ HR =« , F 1 , 17X 1 1 H-*/ i X, X, iH» 

f/I' X, : H>-, r'^y ,■»} -'L T FF‘-S<:uFE (ATM) =« ,F8o1 ,I 7X,1H*/^ ?'■ X, 1H» 

6 / 1 ' Xt.>'.-f 5 7 y ^ ■, 7 X( ' )) 

FT=Fr;-ny+p^'Cy 
C0NV=F- FC/ { F ■■"•:3 + P CY) 

X4I=XAF*C.0’1V+R .CY->r,f;.PIHR/FT 

xpi=xpR^'-cnxv 

yNI = XMF^?’CC.\V 

■'.HI=XHGP.*!i-F ic Y/FT 

XHYI =::HYrv^-RECY/FT-R cy^cffinr/ft 
A = FriLO 

FT=ABl+h..;CY«;n‘'./A 
N=(MWF‘>2c«-)'KF + a.', ^''■^iF 

PRINT?>, (PT( n ,I=7,T- ), {FW( I) ,I = i,IR) 

PRINTftr 

22 F;JRMaT{ 1H1/1' Xt^JT.'y Ac, 9 (P.'tAo) /) 

62 FORMAT ( q-X, 6 i ^ C , ^QBJFCT I VF FUNC«) 

INC=3t. ' 

IXC1=I!1C41 

INC2=IHC/’ 

XN=XNI 
XP=XPT 
XH=XhT 
XA=XAI 
DO 3! 1=1,2 

A' -CPA ( I ) «X A + CPP( I )*-::P + CP;-H I ) ■'•■XN + CPH ( I ) *XH+CPHYn I )#XHYI 
A1=CPA(I)-CF4{ I )+:?< *CPH( I ) 

A2=CPP ( I )-’CFNl{ n 

A3 =C PH Y { I ) ■»» N /3 < ' C P P f T ) - C N- ■ 3 ® ) 4-C PH ( I ) /3« 

A4=CPHY(I )«N/3«'-CPN{ I)-CPH{ I)«fv/3o 
AD(I,l)=A. 

A0(I,2)=A1 
A0( 1,3 )=A' 

AD(I,X.)=A^ 

ADtI,5>=AA 
31 CONTINUE 

DO 25 ITHA=1 ,ITP3 

X4=ro 

X2=ra 

X3=j« 

X4=Mo 
REWIND i. 

DO 5 IJ=1, IR 
P=PP (I J) 

TI=TT( IJ) 

H=WC(IJ) /FLOAT (INC) 

AACIJ, 1 ) = WC( IJ )^-FCrO/ABl 

AA(I J, '!)=TI 

CALL SIMULA 

AA(I J, r ) = TI 

AA(I J,A) = AFT^{ a'.+XA) 

AAIIJ, F ) = AFT^K /:MI‘{X^+X?+X4 )) 

AAd J,6)=AFT^{ X?+X2-:'3) 

A A ( I J , 7 ) = /^ F T ’’* ( U 1 + 2 •s X 2““ ( A” 3® ) /3 o** N* X4 / 3® ) 

AA(IJ,P ) = AFT4( XHYI + (X3 + X4)#A/.?* ) 
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23 


COHiriU"" 

IF{ .'C'-. .U " 

B'L=r:M 

r^ 2 =;M+.; 3 l + ^-I 
B 3 =o 

, -.‘c + 

) 

:+a)/A 

‘> 0 ^ e J' i "* 

V 1 4 f ( T 

( .. X,I 3 f - 


i;i = B 3 ‘* ’ 
N ) / '. ■ 
CG2=(' 
IF(IT'- 
P o I ’'] T ' 
FCFM'X 


P, 


; + /;3«' (33* {N* '3. 


)/ 3 «*“{ 13 »- N)/ 15 o )+X^*(B 3 *N/ 3 o‘-{ 


F ' 1 5 L 5 




A 3 




A I, 


'i £. 


ij 

a 


PRPlT- ,<?. »CiJ 3 
fdrmaT(:. •' f "'F 
GCCl )= ..+'^-’‘3 V .^.:2 
D 0 ( ?, ) = ' 3 '*C 2 
noC'. )=c::'>( i -o*’ n /: 
nD(<i-) = c :*'■( j n / - , , , 

33(1 ) = C 3 >- ;■.+ ('->/.+*%. !)• ’ 

2 C(':)=nr'i )«(3 +- 0.1 

PQ (■;>, ) tiOfj ( ')•!({ Bi + Xi ) 

r)0(/^) = r'D('.O« (3* +■•'! '> 

I F { A 35 ( C '■ ■ C C 3 ) V 3T o ..R - 5 
I F ( A 3 S ( X.,.' C 1 )o 5 

iF(iTFci:w-= ) '30 TC 
IF ( . Cc HP ’ ' GD TO iA 

ni\/=niv: 

CALL nPTKT-raC'M 

IF( ICBcF&ai ) GO TO *1 ■ 

D 2 =C 0 F, 

Di=Xi 
DWC = ' e 

DO A 3 IJ=liIF 

DWC=DWC + Wr'!( I J 5 ^ 

IF(nWCcLTc «) GO TC 

ICA =1 

DIV=DIVA 

PRINT? , 

form at ( 1 ' >' t ■'*F gun C + V c* ) 

;i ^3 ^ = 4 B i 
A.BC=ABC + '^a 

n f j 44 ] J - *t. f I f V ^ 

P J p, 1 4.p -■ C Y S A B L / F F ; D 

GO TO 

DU 1. 6 I J = L » I F 
nu lA 1 = .. 

I C 4 = 

GO TO 

PR IH T 8 

pijp.Mi\T ( 1 ' ' X' i*-CU:'!„T 

ITf= 

Div=ni 

PO A'i = 

TT (I J) ^ ^ 


33*C2* (N’ 3o ) / 3® 
(N* 3o )/3«> 


GO TC 
GO TC 


VIOLATED* ) 


i, 5 o H ) 



HC( IJ) =^■^ ■- ( IJ, „ 5-AR„/F- -D 
CC IJ = ... .!!> 

4e;r)(ij)=,., 

TIC=TT ( I J ) +1 J'i ( TJ ) .^apc/CI V 
IF(Tir>.Gi'.. 1 or TO '.'.i 

IFCTICcLTc .2 ; •> ) GC Tfj ,.7 
C-n T 3 '■ 

/iBD( IJ 5 c- ( IJ ) )^OIV/ {TIN(IJ)*AeC) 

AC2{IJ)= 

GC TG : ■ 

ABO( IJ) = ! ! 2£- tt( I J) )-»OlV/ {TIN(IJ)*AeC) 

AC7( IJ ) = . G i>'o 
CONTI '4 U 
IC]. = 1 
no = 

IF(ACP( J JlcLTc, A3D{ IC- ) ) ‘IC?. =IJ 

CGOirsin 

UO IC . ) 

wC .' = If 

IT’ = 

!iC.,.= o 

DO 0 9 IJ=.l,IF 

TT(IJ) =TT{ IJ J + TI'K IJ )«^\B3^AEC/DIV 
WC(IJ) =WG( IJ ) + wJ'i{ IJ )-*AB“/D IV 
WC.. = wC: +w I 3 { I J )*AP //Cl V 
I F (WC{ IJ ) oGT^' o ) GG TO 09 
IFCAC) .GT. nC{ IJ) ) IT.=IJ 
COHTIMt- 

IF(AR'''tLTo .o ) TT( IC* )=ACr.( ICl) 

wci=wc.' -^-F^ :c/Ar 

WC ?.=■•* wC„, 

IF (WClcLT« :R ) GO TC 
IF(ITlo'-'w’ ) GO TC 2c 
AB<i=-^r, { IT i )/ (WPH ITl)*A32/niV) 
no 46 1J=I,IP. 

TT(IJ) =Tr ( IJ )*’TIN( IJ )-»AB3*AE43A3C/DIV 
WCCIJ) =WC { JJ )".'rt!Nnj)-*A8?»AE4/DIV 
IT ^=2, 

CGMTIMU 

PRINT''.' 

FORMAT! )!,«ITFRA.TI0‘I LF.SS FEED PROGRAMME AGAIN*//) 

GO TO 
P K I M T C 'I 

FOIMAT {//•::, ’.SG THAN TWO TEPPERATUTE ARE AWAY FROM BOUNDARY*//) 
GO TO :,/ 

CQNTIMU:; 
no 9 j j = ,,, IR 

DO 9 1=;. ,c 
AAdJ, n=AAMI Jfl) 

XI=D.. 

C 0 ~ ; 

GO TO i/ 

PRINT.? i. 

FORMAT! ..Mi, OPTIMIZED RiSULTS ARE*) 

PR INTI''- 
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r 




FowM -.1 ( ////: :;icT.]R v^t cata t in r t out r aron'at naphti- 

l PAP4FFI ! H i FYDF. 1C*//) 

OH ' I = . 5 1 p 

PRIflTl ,fAM I, J),J = ^,r:) 

CORHAT F'.Pt ..»r3oA,F9ol,F7«l,F8olTF9,I,Fl, oltFS^l/) 

i.2 COilTIHL- .■ 

?F I IT* .see,, 

-F FOF '-tM ( //. 'f'C .OT.'Sr OF ARCK.ATICS I.N C5+ AT CUTL£T = «t F6o3// ) 

STOP 

IBFTC SUB' 

supROirr I s ihula 


c 

C THIS SUlSUTI’l'- SIHULATO TH£ .^FACTOR AND PUT THF DERIVATIVES OF HOCFL 

C "'Tn-JS O'' UNIT 


CMMMS /BLS -iS'S 

COYHOM /5L.S : » .::.f -s =‘V,*:F tA,;',AtXN,XP,XH, IC,CP,Cl,C2,C3,C4,T 

COMMU'I /OL'S A,.. JrASAy-tPj FT 

COVMCM /.IL S D 

CTPMO'i /BL-S Cru.H 

COPHUM /'3L?/ IJ,FU(:S;’) 

C'JPMON /BLT/ T'lr, p,c '.jINCT 
COPMC'M /8LS/ ABC 

Cn^MOH /PLIC/ !v 

COPMON /3L.S/ TI,XAI,X?I,XNI,:THI,aHYI,A^. 

CaVMQN /i3L'^'?/ AQ(2,0) 

DIPONS lU^ D ( -'i, S Ni ) , (P ) 

REAL H,K{.S3 ) 

Ai; {T)=. ;:p ■•^<'• 7 i «./t) 

A2 H T ) = .■ „ P { A •■:.o 1 A • 4 5o /T ) 

Ai2{ T) =■■■', ,P( 3To9'-.5S6 i o /T) 

A2P{ T) =’U-S{ : , c /T'7S2) 

A,ii,( T) =- S( , /T) 

Fi(P)= j'f ?■'■*{ " *TA^XH*--3’ 3/ (;Z2*A**2 ) ) / ( FT*A ) 

F2(P)= 'S^M'^-.S''!‘*;'.h/A--TP/F4 )/(FT*A} 
f-' {\) = :':r* iP/ if j^\) 

F4(4 )= '3-pl\/ ( FT-^^.A ) 

F3 ( C U = { -Si « '91 e ■: , « +19.'S a *-C 2+ ( N-3o ) *8if-fU *C3 + N#223r-SCe *C4/3. J/ CP 
I C =i 

DO 1 I 1= , UC 

T=TI 

xic=:';:. 

/ 2 r=e:i 

X3C=)C’ 

X4C=X4 . 

DO 2 I 

A = S ''IC-ZZC'-XAC 
XA=XAi+>;:c 
xp=-;pi+:- '0700 

■/.h=:;hi +: c ^)'ic- :eoc- (N-i-a )«-x3C/:i®-N/3o»x4c 
YH=XMI‘ (i:,iC + K2C + X4C) 

-:=A 3 .‘ (T) : ■ 

■USAP) (T ) 

> (T) 


EQUA 


I 

f 

i; 
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(T) 

= (7) 

Cj=F..(P} 

C 7 = FP{f") 

C 7 =F 

C4 = F-f{ '4 ) 

=A 0 ( .) -T+A.l{ ) 

i:..= 4 C{ •: , .:)''T+.\';( : ) 
fi 7 D ( '■ t ) 4 1 + a 1 { ) 


A.''‘ = \D( 

? ' ' 

(v;7 

+An(;: 

A4 = 4 

i){ 

, f , 


+ 4 )(/ 

CP =A 

+ 

^ L * 

' „c 

■*"4 1.'’^ 

JF( I 

e 

f ' ' 

7. 3 , 

! C 

-LL D 

K( It 

.-) 

=C v. 



K( I, 





K( It 

3) 

™ C 3-' 



lU I, 

n ) 

= C4 



K( I, 

1 ) 

= F5 1 

; c 1 

) 

IF ( I 

0 ' 

Cc -1 ] 

1 G 

■! T'-) 


■.],C=’7 +h'..K ( Isl)/C(I) 

'2C= P +H«'-K( I, 2 ) /C( I) 

/.4C. = :(:.+tlii'K{I,7)/C(I) 

y, 4 C=:<^’-+if-K{ 1 , 4 ) /CU ) 

T=TI+H'^Kn ,5 )/C( I ) 
cri'iT i-iii •: 
nn 3 i=.;, ,>5 

PHK I) = (K{: , I ) +K( 2, 1 )^^2,+KO, I )*2c+K(4, 1 ) )/6e 
Xl=x;i+H^'PFU ■ 5 
X2=X2 + H*-FHI ( ’) 

Xi=':''.+ H+P|il ( .) 

X4=)<:-t‘+H*PH J {■'''■) 

TI=T I + h*PFT { j ) 

CONTMUi: 

'<a=::ai+::^ 

>;n= 7/'JI- ! .::. + X2+--n 

XH=XHI +' )*X3/3«.4s»X4/3, 

XP=XPI+.7.'"'.:2 
T = TI 

22=A2? (T) 

= (T) 

^^■=A2? (T) 

( T ) 

C?.=F ■.{ P ) 

C?,=F.'( P ) 
c:- =F ;■( \ ) 

C4 “F 4 ( ^- ) 

=ADl 1, X )«T+AD(2,i) 

A.i. = AC( )->T + A0(2»2) 

A'd=AC( 1, '.)'«T+An(2».2) : : 

A J=A 0 ( . J U »n+AD{ 2 * 4 ) 

A4=AD( : » A )*T+AD { 2 j *') 

CP = A ^ +A2«r.2 + A3'»X3 + A4«-X4 

CP=CP+A'' 




FU(IJ,M=C. 

FU(IJ,,';)=C2 

FUdJ,? )=C3 

FU(IJ,^',)=C- 

FU{ IJ, • ) = F:; (C,. ) /ABC 

CALL D^FI 

WRIT.; C) ( ( (D{ I, J, JM) , l=i,5) , J=i,5), JM=1, INCl ) 

R'lTlJB/' 

= ND 

I8FTC SUB2 

CUBPHUTIrJ-' C RI 
C 

C THIS "Ur.-'nLTII CALCLLAT.'S CFRIVATIVES OF FCDFL ECUATICNS 

C 

CO-WMCM /eL2/ ‘ .'J, :<','’5,«tDD,e,FtC,ICtCP,Cl.C2,C3,C4,T 

CO^^MCN /BLm/ '■) 

Cn'^MOH /;UV A '\2 '■ FT 

COMHQM /«LC/ 4’3C 
CnwMOM /BL.,'’/ '! 
ccyMU'-i /tl:7/ 

DIMFMSI'VI B { ) ,G {5 ) 

P. L AL A! 

ic=ic+:. 

B:^=rT» p/ { FT-* ) 

3 4 =p »» 3 s-C ^ 2 / ( ;V>-» r. "• A ) 

0( I, IC)=Bl->{- A-*(C-H9« *DD )“12»*C*DD ) ) 

D( It lit IC) =Bl-*('- A + 8' 34*{““3e«•A*CD■^4«*C*CD) ) 

D( it 3t IC) = S1#B4-»^CD»( N- ’o )«A 

n( It 4t IC ) =ei^‘( (-,N-»DD^^A■^4*»D0«■C) ) 

D{ It C, IC) =B!->A^B#.V''7 3i „ /T»<-Z+P««4«DD*C*«3»EXP ( 1I295«/T“22*94 ) #11295 
:i9o/( A»*4^'T-»#"-*-FT) 

8I=„' »P/ ( FT-fl-A-»#L ) 

84=P/A 

n{2t 1 1 IC) =31*( 'C+C«*B«A•'6o*B•»C)■^3a^|•E/H4) 

0 { i t 2f I C ) = B I-**- ( B '■•■■*■(’- A C" 3# 4 + Eo’-B^-C )•”{A■^E) / £4) 

D{ 2, Bt IC ) (\i. -Jo .o-^A/!E43 

D(2t4t IC)=e2->(B4-''d- A-sC-N^-A^B/Bo + ao *8-»C)“E/E4) 

D{2tFfIC)=R2’^P-s-B'*'C-»596' • ^/T^-^I'-CaPI 43 « I “’67 6 jC o /T > «P-*'F.-*'67 / ( T-»*2 

i#FT#A) 

R] =’::5/ { FT-’s A# ) 

DOf It IC ) =• El# 

D{3t2t IC)=S1*(/U':) 
o{ It 3, ic)=’-'B:i»a 
Ddt^t IC )=BI-4-: 

D(3iF,lC)=ei#A-s-c«-&.'3‘ , «/T##2 
0(A, 3, IC)=B1 *('-‘V>Ik.^8) 

D(4t 2t IC)=B1#(*-A-+B) 

D { 4t i't IC ) = - « 

0(4,4, iC)=31#{' A + B) 

D 1 4, 6, I C ) = 3 1 '*'4 ^ B#6 ’’S t, e /T** 2 

DO 1 j =:i.tc : 

X G { J ) - ( * 0 ( < , J 1 1 C ) # 9 i 5' <1 e +0 (2 t J f I C ) »19€<!0« ( N'-’S « )*D(3,J, IC)*3i*A e-*-D( 

■4, J, IC )#N#7433« )#CP 

81=-C^ ^’■9; f '•%+C3*l9i:.,!,o+C3* (N*--3« )-»8lCC'o+C4#N#743:3® 

nr-:, *» ic) = (G(i )-3i#Ai)/cp»*2 
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r( If, ) = ( G ) /cp4>'<^2 

0 { » .i j I' - ) = ( G ( • 3 1 ■ . ) / C 2 
0( 5, 4, iC) = {G (-■■,)■ 2 

DC;, =, JO =G{''C/C?«^ :> 81*(AD(i,l)+AD(l,2)+AC(l,3)+AC{l,4)+AD{l,0 )/CP«*2 
D0 2J = .,^' 

D( J, 1C)=G( J,, , !C l^’OG 

2 c( ", j, ic )=r; ( j, ic )/’. 3 C 

PP.TURM 

cm 

IBFTC SUB3 

SUBRGUTCi L'Vli 
C 

C THIS SUB-^ntTC! ' IH.OnATh THP SYSTEM VARIATIO'^ML EQUATIONS 

C 

Ciwrc.m /&L';/ D 
COMNUM/nL • / LI 
CaMH'jM /JL.S/ CC'),1- 
CUMMUM /•JL3/ I ir , li-c ,T\C2 
R M. L I ( O , L ( ; ) , i< ( A , ) 

CIMEMSM 1 u( i, ,50 
OIMCI-jSION IA(^) 

IC=INC;, 

IA(l)=r 
IA{2)=j 
, IA(3)=’ 

' IA(4) = 2 

DO 8 11=1, INC?. 

DO 7 1 = 1 ,1 
7 K I) =LI{ I ) 

DO 9 IJ = .,,4 
ICl=IC'It(IO 
DO I* 1 = 0 5 
K{ Tj,n= , 

DO 1^: J = ;',5 

C K( IJ ,I j=K( IJ,1 )■ L( J) •*D( J,I, IC .) 

IFdJolQ.^O GO 'C 9 
DO 1.1 1=.,1 

IF(ABS{L( I) JoGTcLp +• 9) GO 10 15 
Cl L(I)=L J( I)+F^K{I.J,I)/C( IJ) 

9 CONTINU.: 

lO 12 

i?. LI (I ) = LI { n+H'» (KCi.,! )+K(2, I )*2 o+2b*KI 3,I ) + K(A-, I ) )/6o 

3 IC = IC-.': 

Rr.TtJPN 

15 PRIMTI • ■ 

i F'lRMAT ( VALU:S ARE OVER FLOWING IN LAMS*//} 

STOP 

ftin 

IBFTC SUB4 : 

SUBROUTl'T. LAHA 

^ I THIS SUSkCITCH is TO INTEGRATE THE .ADJOINT EQUATIONS 

C 

CnvMCH /i?.L4/ D 
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COi^MCM /BL!;/ LI 
CGMM'T'I 

CfJMMOM /3L?:/ I -’C , I NC L , I NC L 
Rf’AL LI(=;) ,L( ) 

Clf^ENS JHrl D ( -t . ) 

DiM- jsn-i I. (. i 

ic=i 

IA{1)= 

) = . 

I A ( 3 ) = 1 
iA(4) = ,: 

DO 8 T I=. , INC2 
CO 7 I=. ,3 

7 L{I)=LI{I) 

DO 9 IJ = ..,4 
ICi=IC + lA( IJ ) 

DO I' 1= 

K{IJ,I )= , 

DO I J = 

i:? K( IJ ,n=iu IJ, I ) + L( .j)'‘D{ I , j, ic: ) 
iFdJoiiCoA) GO TC 9 
DO 11 I = ’ , 3 

IF ( ABS { L( I ) )oGT3^„ ■■+ 9) GO TO 15 

11 L{I)=LI{I)+F*K{I J,I)/C{I J) 

9 CQfJTrjU!:' 

DO 1? 1= .,3 

12 LI{n = Ln I ) + H*(K(. , I )+K{2, I )* .<,+2**K(3,I) + K{4,I) >/6« 

8 IC = IC+2 
RETURN 

15 PRINT! 

1 FORMAT {//..(,« >/ALl ;S aRB OVER FLOWING IN LAMA«//) 

STOP 

END 

IBFTC SUR5 

SUBRCUTIN” CPTI{TT,i:3) 

THIS SUBRCUTI^l-: CALCLLATFS GRADIENT VECTOR FOR DECISION VARIABLES 

CQMMCM /BLI/ X:,,.:2,3:,XA 
COMM CM /rsL‘'i/ D 
COMMON /BLll/ LI 
COMMON /2LS/ C ( *' ) , H 
COMMON /3L7/ I JjFKG,?) 

COMMON /BLS/ I NC , I NC,, , I i4C2 
COMMON /BLIP/ j 

COMMON /8L1../ DD(9),C1 ,C2,E1,B2,B3,CIV,IR,IR1 
COMMON /BLi’/ T!Nd),WlNC5) 

COMMON /BLi4/ wC,A 
DIMENS lUf'l WC ( 3 ) 

DIMENSION TKS) 

CIMLNSIOM 0(D,9,r.i),nH 6,e,3) ,AA(2,3) 

REAL LA{::,5,9) ,LS{ j,:., 5),LA4(3,5,2I ,NU(2),LI(5) ,N 
REWIND 1 
DO 2 I J=l ,IR 

READ (I) { ( (DC I, J,JM), l=! , 5) ,J = 1,5} , JM=l, INCl) 
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H=:«WC{ I J )/FLr)Ar { IhC' ) 

cn :: 

DO J » 5 
4 LI ( J )=? « 

LUI) = :,„ 

CALL LAMS 
DO 3 J=.'. T- 

2 LSdJtl, J) = LI{ J) 

2 COHTIMUr 

ROWING 1 
DO 6 I J = l, IR 

READ (1) ( ( (n( J, J,jv), 1=1, f ) ,j=:l,5) ,jf/ = i,lNCl) 
H=WC{IJ)/FLCAT( IOC ) 

IJK=1 

IFdJ^' Q,! ) IJ.<=;= 
no 7 I=IJK,5 
DO 8 J=:i 

8 LI(J ) = ■ o 
L I ( I ) = 3. o 
CALL lama 
00 7 J 

7 LA(I J, I, J)=LI( J) 

6 CONTIMUF 

LAA( IF f I, I ) = ;c 

LAAC IR,:,/) = ?o/A+(B' *P< /A^^2 

LAA( IR,2,1 )=■• o 
LAA(IR,2,2)=* (31 + X;, 

LAA{ IR,3,2)= e 

LAA( IR , :.'t 2) =' { B1 +X_)' (“do+O/i ?e + { )»B3/3o ) /A*»2 

LAA( IR t4, i )=' a 

LAA( IR >4, 2 ) =•'( Bl + Xl ) # 'o +N#B3/39 )/A«*2 

LAA{ IR,r,i)= ® 

LAA( IR , 5, 2) = > 

CO 22 JI=1,IRX 
IJ=IR+1-JI 
DO 9 1=1,5 

sc=f:u 

SD=0® 

DO U. J=l.,4 

SC=SC+LAA( 1J,J,1)«LA(IJ,I,J) 

10 SD=SD+LAA(IJfJ»D)*LA(IJ,I,J) 

LAA( IJ'-ltljl )=SC 

9 LAA(IJ*‘1iI»is)~.>D 

LAAC IJT5,i5=LAA( IJ* :.t3f D 
LAAC IJ,5,2) = LAA{IJ- „,5,2) 

22 CaMTINU'* 

SC=i’o 

SD=&8 

DO 13. J=lt4 

SC=SC + LAA ( I » I ) -^LAC I »5 » J ) 

11 SD=SD + L.AA( 1,J, 2)#LA(I ,0, J) 

LAACi».>Tl)-SC 

LAAC 1, 5,2) = SD 
DO 12 IJ = lf IR 

S D “2 ® ■ ' ■ ' ■ ' 



J-C= 

C 1) ' I = I f 

^■C=SC+F-J( JJ,I ) -LAA (IJ,I ,1 ) 

i:, sn=Fn+ru(ij,i)-M.AA(iJ,i,2) 
nidj,;:,: )=-^sc 
DI ( I Jf ',<■)=” 

Cl (iJj *. f d 

DI (J Jt =L4A ( TJ, It. ) 
ni (IJ, =l ( IJ, ) 

12 [)I (TJ, T'')=f 0 

I c :! = I R 
TC2=': 
ir,2=,i 

f'J no 1/: 1=.,/; 

DU I, A J = ,.. 

14 D I ( :> , I , J ) = L S ( , I T 5 ) * 3 1 ( , 2 , J ) + F U ( I , n * C I ( I , 6 , J ) 

Oil IS 1 0=1, 1 R 

DC iS I=,.,A 
SC=' c 
SD=!> 

St 

CO 17 J=l,« 

SF=LS(1JtIiJ) 

SC=SC + DI( IJtJi iJ»SF 
SO=SD+OI( IJ , J, D^SF 
5i:=S&+DI{ IJtJ, lO^SF 
17 COML-lljr 

SF=FU{IJtI) 

DniJ+ > , I ,1 )=SC + DI (IJ,6,1)*SF 
DI(IJ+i,I,2)=SC+DI(IJ,6T2)*SF 
DU IJ + lt i» J )=S,r + DI (IJt6t2. )*SF 

15 COMTIMUf: 

AAdfl )=nU IR + It It i ) 

£A(1 T/;)=OU IP + 1t It 2) 

AA(1 T? )=''CI { IR+ItI t'U 
SC — . o 

S F ® 

no If 1=1 T 4 

sc=sc+ni( IR + I, I,.-) -^DCd ) 

Sn=SC+CI (IR + i T ItO'^CCd ) 

se=s:'+ni hr + It I f 35*DDd ) 
lo CGMTIMUL 

AA{?,:)=sc 
AA(2,2)=SC 
AA (?.,?’)=• SC 

[)T = AA< It J.)«AA( It?.)- A A ( I r 2 ) * M ( 2 t 1 5 
IF (Die V ) GO TO IS 

NU(l)=(AA(it3)^AA{2,2)“AA{2,3)*&A(l,2) )/DT 
HU(?.) = fAA(lTl)'“'AA{2T2)‘”AA(2*l}*AA(lT3J)/DT 

D*"! 2 4 I J = I T I R 

DTK-=DI { IJ t5 »1) . )+DI { IJt 5t2 )*NU( 2) + DI ( IJt5,3 ) 

D'«C = ni ( IJTSt’ )*OU(I) +DI {IJt6t2)*NU(2)+DI ( IJt6t3) 
TIN(IJ)=CTS 
WIM( IJ )=D'aC 
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24 COMTIMUr 

I FdC?. ) R TUR , 

DO 6 I J = ;; 

IF(TT( IJ), C,>‘' ty) GO TC 61 
IF(TT( IJ)..:S3 o 5 GO TO 62 
GO TO <■. 

6T IF(TIf!(lJ).L:;c ) GO TO 6 
GO TO or. 

62 IF(TIM( IJ)oG’;o o) GO TO £r 

63 

Q!(IJ, 

DI(lJ,S.d='„ 

ICl=ICd 
6f COMTIHUL 

IF{IC:;d*C«IR) o: TtG'l 
IFdCX.G'o?) GO TO 
IC3=-, 

RETURN 

64 IC?= 

GO TO 66 

18 PRINT,?' ^ 

?J FORM&TI/ZlXri^ il 'Tr ‘^MINOr^T IC lizPO*//) 

STOP 

FND 

IBFTC SUB6 

SUBROUTINS CPEV 


SAME SUBP.QUTINf. 


AS GIVEN IN SIMULATION PROGRAMME 


ENTRY 
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JOB 
IB JOB 
I3FTC 
C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

e 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CGG:-., ■' 

P-JOGP. A r-CP 


9 f ? AG 


P jNAME -^SHOK KUFAR SHARMA 


,'T !lv3Z^TlCil CF OPERATING COST 


IR IS .lUilG P 


7! F 


rADTMv T- vn- r"'. A^^^D^HIGHI-R HYDROCARBONS IN RFCYCLB GAS 

P i-'-S' CV ‘"Ar. ^nL-s CF i»PriTHS ffcd wo rccycl. gas pw hour 

pp'tc zt. -'F-.p, -,C I 4 L T STR'-iAM TC VARICUS RFACTORS 

TT AND v-^c’pR : T-!!p': ■■ TU- S CF INLST STREAKS AKC WEIGHT OF CATALYST IN 

in pi WHICH A R=ACTCR IS DrVIOUO FOR INTfiGRAT ION ; 

\a\ -|F vj G CF 'APHTHA WHICH APS US7 0 AS BASIS CF CALCUL AT 

r- Fcr. CEJcCTIVn PUNCTICN 

ifiu- r"' J'HA’G'f 'T uring iterations at suitable places 

THOUGH TH ■ ^ V 'i ST/ P'^I jG VALfS HAVF B"EN FOUND TO 3 r 

AOVI SF C 1C T^Y OTH T' START! HG VALUES FCR THESE FACTORS TO KINIHUt TOTAL 

numbe r OF ‘ ^?,/p .V prp F=fd TO FIRST REACTOR OF THE SERIES 

DIV IS 1 F'iCTOK jr D CI.V- ^/G^UU^y: CF CHaNG-S IN DLCISICN VARIABLES 

fTR^H. Ja'I^UM Jr-UTlSNrFn^WHICH THE PROGRAMME IS ALLOWED 

FRi''aND^HPJ aRc AMFUJTS OF errors ALLOWED +N CCNST-A + NTS ONE AND TWO ON 

'^’^oLrClIV^^^^ function is less THAN FR PROGRAMME IS 
TeRMINATHO*' AUTOMAT I^•-LLY ANC THE SOLUTION VECTCR AT THIS PLACE IS PRINTED 

IT ?rroJlSEo"THrT"M.MY INITIAL SQLUTIGNS BE TRIED TO ENSURE THE VALIDITY 
OF THr R .SULTS 

mole FRACTIC'IS OF f-’CYCLH GAS COMPOflENTS STARTING FROM METHANE TO HEXANE 

AND HYCROGCN I I -'FO 4 FORMAT „ « nurucMrc f r rn rcn AMD 

WEIGHT pr.r< C5 1T OF PAF4FFrNSlC6 TC C9)? 'N4 PHTHi^Hw^ ( C6 TO C 

and' AROMATICSC TO CS) IN lr F8®4 FORMAT 

T?Sp!:^UUp'':'v^'Jtmh,.1^'‘iGHT of catalyst VECTOR AND PRESSURE VECTOR leFlLct ) 
Mo'lF FRACTICNS OF Go A'lO HIGHiP COMPONENTS IN RECYCLt. GAS (l.F8«4) 

FEED AND RECY t''F3 


D AS A STARTING VECTOR FCR OPTIMIZATION 


ju^-cc: i-i- wC VALU''S AR*- L- r, - - - 

jMfof VALUES A,LL THo CONSTRAINTS SHOULD BE SATISFIcD 

COMMON /BL]/ XG.»N2jXi»X4 
COMMON / BL 3/ A ^ j A2 » A"* » A4 1 P * FT 
COMMON /BLc/ C (O jH 
COMMON /3LT/ IJjFU(5»5) 

COMMON /3L8/ I HC , I NC 1 » I NC2 
COMMON /8LS/ ABC 
^^MM0^4 /BLlC/ 74 

COMMON /3Lil/ CD(4)»CX»C2»Bl»B2»B3»DIVTlRtIRl 
COMMON /BL’2/ TI »Xi. I f XP 1 1 XN I tXHI » XHYI t A* 
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!*■! /3L '1/ tt'K -) ,WIN(5} 

/2L.A/ 

CC‘'7C.t /3L; = / C“>A{ .') ,CPN(4) ,CPP{2) ,CPH(2)»CPHY(2) ,CPHYI{2) 

/9L:i/ PP,>rjF,XAF,X!-Yi<,XH2R,MWF,FWR6 

c /3L'7/ sn(^,n 

COYMOH '/3L?'-V /.JUT 

COWMC;'! //LCl/ 'IX 

CI''!''''IS n ! FF(" j , .-^'C ,TT(5) 

D1.^^--XSIG/ IM", ) , ^ , FT (5) ,P1<(5) 

ni p"x:^ 1 / i a !P c ' , ) , tout ( £n ( 5 ,6 ) t xout { 5t 6 ) 

■jp'rxciin .10 'K ) a " a :) 

P;:AL ■ 

DATA , -.H T£ t JH T3 ,SH T4 ,6H T5 ,6H Wi , 

1 6H W2 W • 5 .'H fl*-i f SH W3 / 

CALL CP V 
4B.l = U» 
a32= '• .. 

A8C=.fP 

T' =12. • « 

Fft = '» c 
'•Rl=c' ' 

"•:R2 = c''' i 

ITR3=5 

DIV=®.?.S 
01 VI =D IV 
DPAO A: .IP 

READ 61,(71(1), WC(I),PP(I),I=,IR) 

61 FORMAT (8F:1 «: ) 

60 FORMAT (It) 

IRl=IR*'i 
SC 1= lo L +3 2 
C(l)=2^ 

C ( 2 ) = / o 
C{3)=1«. 

RrAD61 ,F :i-;C ,PLCY 
RCADol , C6FINR 
AFT=Fr': D + POCY 
A = FEED 

DO 1 F I = , I R 
IS WC{ I ) = WC( I )«An.../A 

PRINTS 
PRINT2 

3 FORMAT (l.Hi) 

2 FURMAT( ////!'■ X,"'2{ iH« ) /IC >',,1H^‘,T0X,1H*/10X,1H*,30X,*FEED DATA*, 3! 
IX, IH*/:' ■',1F>,7' X,72{ lH*)/lC?X,lH*,7f X , 1H*/1L X, IH*, 7t X,1H*) 

PR INTI 7, F~SC iRFCY , MWFjMVvRF 

17 FORM AT (£,.■ X , 1 , 2< X , ■♦('APHTHA F'.ED LB MOLE/HR =*, F8e 1 , 17X , IH 

1*, IH*, Z1X,*RFCYCL1 GAS LB MOLE/HR =*, F3o 1, 17X, 1H»/ It.X, 

21H*,7 /, .F*/V. ::,:h^,1?X,«L8 molecular wT of NAPHTHA =», F8« 1, i7X, 1 H 
3*/i:.X,lH*,7^'X, iH*/iX X,IH*,12X,*L8 MOLECULAR WT OF RECYCLE GAS =«,F 
48ol, nx, lh*/n/',lH-®,7uX ,1H* ) 

CD( 1 )=::AF*FS'‘0+RECY*C6PINR 
CD(2) = X'JF«FcnD 
00(3 )=APF*F=fO 
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*^'CY-RiCY*C 6 PIHR 


, o 


62 

22 


'\i It.' 


) 


cc('';-) = ':hy 

«r CY 

F.'trfUT ( 

H*/' ’/v-P', LB KCL‘/HP =*'F3altl'7ATiH'S'/l-X»iH*»7*A»lH*/?u 

’■ ,«!'Y0^-1C^BB0MS LB '‘^'BL../HR =* f F8® 1 1 17X) aH*/ li X » IH* t 7‘ X» iH 
■ h-/ , '5. -' j-.-HYBR CG*^ '1 LB MCLE;/HR =«• , F3o 1 » 1 7X 1 1H#/1J-X » i H«7i X » IH'**' 
, iLYT PR-SSUFE (ATM) =■» , F8 « 1 » l^X , 1H»/ 1 ' X, aH* 1 7* X , 1 H« 


V’j IH’‘»"'*X»->"''.RDMATICS lb MCLE/HR =*TFSol,i7X,lH«/i<''X,iH»,7 
•'sGH-"-, :■ , /-’iApHTHtNES LB MOLE/HR =«, FB« 1 » 17X1H«/ 1'_X1H*7(. XI 


. 77 ( IH^ ) ) 
CY) 

'CoPIMR/FT 


> ■: 

’■’' ? / i j 
‘3 / -I . . » J. F ' 

: / ’ ri-;.' J ”• )i » ' H-'* / 

FT = F.' ' G .:y 
CnFV=P n/jF'.'-G+R,; 

.‘U, 1= .;.''.F-Ci!f' V+R GCY^ 

XPI = XPF^^'CCF'« 

'f.ji = vt-iFv-cr: '-v 
::Hi = YHrp->p> cy/f^ 

XHYI =XFYP-'.P”CY/rT-R = CY>C6PUR/FT 
A = Ff- ,1 0 

FT-AB'+PCY^ASt./- 
f.; = ( ’(((/ F ■ . ' 3 •* . K F + ^ ' f. ■*• A A F ) / 6' o 
PRIMT2G, (PT(I) ,T=] ,IP)» (PW(I) ,I=ltIR) 
PR IMT<! 


. f c < 

/I 


*• C i t b ' 


:,»C?*,AX,*0BJ£CTI VE FUNC») 


Xt-*t-IT ;P #tAh,9(,.. tA6)/) 


FOP MAT { P' , 

COPMAT ( ..F, 

XN=X''iI 
.bP=;':Pi 
XH=XHI 
XA=XAI 
DO 31 I = j„»A 

A =CPA( I )'K ’<A+CPP { I )^-YP+CPN( D-tsXN + CPHC I )*XH+CPHYI( I )-»XHYI 
A1=CPA{I)- CF?1( n+:’..*CPH(I) 

A2=CPP ( I )-CF'l( I )■ CFHd) 

A?=CPHY(I)-»N/7j' CFP(I )“ t N*'3e ) *C?H ( I ) /3*, 

A4=CPHY( 1 )^N/3c-CPN( I )»CPH{ I)»-N/3» 

AD( I 5 = ^' 

AD( I t2 ) = b.. 

AD(I,?)=A’ 

^0(1,4)= A? 

AD(I,F)=A*'i 

CaOTI’IUF 

IT!-= 

IMC=5 

IKC1=INC + 1 
IOC 2= 

CD ?.L ITR'4=1 ,11'' 3 


Y —ii ‘ 

’b *, — ^ - 4 ? 

*.? 7 ;, 

Y — c 

Xi'= I 
Rr'ft'IMO 1 
DC *5 IJ= , IR 
Ml AIR { IJ ,?.)=)< 1 
aIOR ( I J , 2 ) =X 7 
::iNR (I j, 2) =:<3 

XIOF (T J,4)=Y-V 
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P = PP ( T J ) 

TI=TT( IJ ) 

H = ’WC( i J) /FLC* Tf INC ) 

; (I Jt . { ] J D/A3,, 

J, )=Ti 
CALL Cir'-ULA 
T']UT!IJ)=TT 

laCIj,/' +::\) 

i A { J. J? ) =A FT- ( ’ !I‘ . +X.: + XA)) 

' JJj ' ) = '-PT ■ ? '3 + '' ‘ : 5 

A ( I Jt ■' ) = '• FT' ( h-'lU )* X3/3o*'N'*X4/3o ) 

AA ( IJt • ) =AFT-s' { ;i-Y I + ( .■~ + /.4 o ) 

5 CO?iT jr-'i,]., 

TF( I TS"-" c : ) C. =/... 

31=;'. U 

b:.='kat i + 

' P I ■ ' 

j>5 + .'-v-R + ,-!*■( B ^ i*’ ( N' 3 c ) /3e •“' ( 1 5©*’ K ) /!'' » 5 +X4'**’ ( B,5*K/3 o'” 5 iPo N ) 

ca< = (F" +:/ )/- 

IFC ITRAc ..Go 1 ) C.,=CQ2 

PR IMTC ;,ITF.4t?TT(IJ),IJ = i,TP),(AA(IJ»l)»IJ=l»IB) 

23 FQRM/T(1 X , I " , 1 F8..-1) 

IF(AP' (C.- X: 5; GTc .'fi' ) GC TO 6 
IF(ABS(C:.‘'CC:') -CTc ,R2) GO TC > 

DIV=DIV1 

lF(TT^;o''Uc 1) on TO . 

rF<IC4t^ : Qc ; ) on to 4^ 

CALL OBJL (T 1 , T' KJT , :.(I NR , ABI. , AB2tTL ,SC) 

SC=SC»-A3'''»'rAFT 
PRINTS-' ,.:it C'1 it SC 
4\ FORMAT ( 9!;:.;'! , 2 FI c4,>:„*'-o4) 

IF(SCcGTo SCI ) CJ TC 47 
I F ( ( SC ' SC ) a LT e !R ) GC TO 2* 

GO TO a;: 

47 A8'l = AB„ + o A 

ABC=ABC + 2c 
DIV=CJV. 

IC4=I 
PR INTT 

1 FORMAT (1 ' At FflUiNC +VE* ) 

DQ 4A IJ = lt IP 
TT ( I J) =A 5.A( IJt ') 

44 WCdJ) =AAA( IJtl )^A91/F7cD 

FT=ABj +r 'CY^ ABi/FcLD , , 

GO TO ::4 
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^7 SC1=SC 

46 CALL 13J' C{T",'^"'UT,;;T;19,aB1 ) 

IC4= 

r,D(i ) = :.,o+?-.^ jc-c j 
CD^ 2 )=--3'-^-*C2 

cc(3)=c;-; M >fC 2 ->( N- ?.» ) /3® 

DD('^) = C?--{ ffj' 'Be } /3a 

DD(’< ) = C2- 4+ (T\+ ' , ■'DC „)•■_») 

CD<2) = C:r(: ) 

CD(3 ) = ")[)(',: .) 

CD{4) = nD( / )« (B‘ +.: .5 
if(Itr.;'.o 'CcI'k-.) g.j tc c* 

CALL nPTKTl? ^C‘ ) 

IF{TCH^:w. ) GQ T- ,r- 
D.2=C0? 

Dl = Xi 

DO Ih IJ = : , r- 
DC 16 1 = :. , 3 

16 AAA( IJ,I )=AA (T J, I ) 

GO TC 

6 PRINTS 

8 FORMAT!] . X ,.-CQ'!ST y:CLi',TrD«) 

IT5=‘. 

CIV=CTV«2o 
CO 46 IJ=1,IR 
TT(IJ) =AAA{ I J, 2) 

45 WC { IJ ) =A.\ A ( I J , ) * A B . / F ~ : ■ D 

41 DO 35 IJ=1,IR 

A8D{IJ)=a, 

TIC=TT ( IJ )+Tri ( IJ ) ^A'^C/DI V 
IFITICeGTalA i; c ) GO TC 36 
IF(TICoLTa 1262:, ) GC 
GO TO r3 5 

36 ABD! IJ ) = ( i4€' c' TT ( IJ ) ) *D I V/ (T IN ( IJ ) *ABC } 
AC2( IJ )=.4fci, ^ 

GO TO ,'15 

37 ABD(IJ) = (126 « "TT { IJ ) ) ^*0 1 V/ {T I N ( IJ ) *ABC 5 
AC2I IJ ) = 1 26.\, 

35 CQNTINIJ*: 

IC1=1 

DO 38 I J = 2, IR 

IF(ABD( I J )« LTb Ajni ICT ) ) IC1=IJ 

38 CQMTINU'- 
A83=ABD{|C1) 

ACl=i'4 

IT1= 

DO 3 9 IJ = .,:.,IR 

TT(I J) =TT( IJ ) + TIN( I J)*A3:WAEC/DIV 
WC(I J)=WC( IJ )+WlN( 1J)#A33/DIV 
IF{WC( lJ)«GTt' » ) GC TC 39 
IFCACiaLTcWCd J) ) GO TO 39 
IT1=IJ 
AC1=WC(IJ) 

39 CONTLNUi- 
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!F( -,n''.,L7wo ) TT{ ICl )=4C2(IC1 ) 

IFdT'f ) GO TO 26 

-6=’^= ,^r( IT. 16 (WI 1( IT. )*A33/DIVj 

DO ij=:i , IP 

TT ( IJ) H ID )-TI N{ IJ H4 33i«-AE4>43C/DI V 
.vC(IJ)=>;C( iJl-i'ilH IJ )^t4B3*AE4/DIV 
ITI = , 

CuGT i 'iU ■ ■ 

LIT' ,1 

FCPM/,T(//-. < ,4' IT.fi; liTIO''! LOSS F‘:DD PROGRAMME AGAIN*//) 

GG TO 
L" PkI’IT' 

..D FOPMATC/// ‘SS THAN Ti^G T'5MPERATt'R£S ARE AWAY FROM BCUNCARY «) 

GO T::i ,. 

w CO-MTnU' 

CO 9 IJ = . , IR 

OC ^ I=',C 

9 AAdJ, n =A:U ( IJ, n 

■<i=0 1 

ca-'=D:' 

GO TO „ 
i PRINT-'' 

21 FORMATC 1H..,:'.,-»;JPTIWI2"D RESULTS ARE*/) 

PRriTV, 

14 FGkMAT(////1 d^REACTGP KT CATA T IN R T CUT R AROf^AT NAPHTH 
1 PARAFFIM H., HYDROC*//) 

DO 12 1 = ., , IP 

PRIMTrdIt(A-'.( l,J),J=Ar3} 

1? FOPNATU X,I> ,F',2e: ,F£d.,F9ol ,F7,1,F8 o),F 9.1,F1 ol*F8al/) 
i?. COflTIMljr 

P R I 'f'l T ‘ ® C 0 ' 

15 format (// i /,*PFRCTNTAGE QF ARCMATICS IN C5+ AT QUTLET=«, F6«3// ) 

STOP 

"ND 

IBFTC SUSi 

SUBROUTI d SIMULA 
C 

C SAM'.' SIJBRGUTIN' AS GIVEN IN PROGRAFME FOR PINIHIZING WEIGHT OF CATALYST 

IB FTC SUB 2 

subp.qutin: C'.’Ri 

SAM’- '^UbnrjlTIN AS GIVEN IN PROGRAMME FOR MINIMIZING WEIGHT GF CATALYST 
IBFTC SUB3 

suBRnuTi.i LAM 

SAME SUBPQLTr-!’; AS GIVEN IN PRCGRAMMiE FOR MINIMIZING WEIGHT OF CATALYST 

IBFTC SU84 

SUBROUTINE lama 

SAHF SUBRCUTINL AS GIVEN IN PROGRAMME FOR MINIMIZING WEIGHT OF CATALYST 
IBFTC SUB3 

SUBROUTIM: CpT!{TT,IC3) 
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COMMON /3L1/ 'Alt Alt A t 
COVMn’'J /3L4/ D 
CC^^MO’i LI 

/?lf:/ CC') ft 
CI^MOM /3L’’/ 

CGN'MCt.' /3L '/ I !C, INC, tlGC'- 
cc‘^MG'1 /FL ■:/ 

CnM*r>M /ULi:/ 3D( ) ,ci,cr:, e:.,B 2,B3,DIV, IRtIRl 
COyMCM /OL:. .:/ TIi( G,ftTN(") 

Cov^nn /liLlA/ >'S.tl 
CQVwirjj /3L,;':/ .'jUT 
COMMON /2L2I/ NTM 
DIML:..iSr'!N '-I'H ,tG) ,ACUT{f-,f-.) 

DIM.-liSin’l -j (N, --I (C,c, 2 ) jaA(2,?) 

DIMC'iS IN 1 ) 

DIME-.'IS inr, TT( " ) 

Ri-AL L - ,2) 

R AL L A f 1 f !: T '"■ ) » L ^ ( t ,5) t LA A { : , 3 ) , NU ( 2 ) r LI ( 5 ) , N 

RCwINP 

DO 2 IJ=' tIP 

RfiAD ( .' ) { { (D ( T , J s Jr* j » 1 =- » £ ) » si=i >5 ) » JM = 1, INCl ) 

H=-WC( IJ ) /FLOAT? INC2) 

DC ?> I=L»''f 
DO 4 J=:,:. 

4 LI{J)= o 
LI{I)=ic 
CALL LAMS 
DO B J 

LSdJt I,J)=LI( J) 

CONTINUr 
REWIND 1 
DO ft IJ = ,\,1P 

READ (1) (({[)(!, J,JM), I =l»£> »J=1»5),JM=1,INC1) 

H = WC(IJ) /FLOAT { INC '^) 

I JK=1 

IF(!J»LQo^) TJK^'-' 

DO 7 I=IJK,£ 

CO 3 J=i.,5 
LI ( J )=• c 
LI (I 

CALL LAMn 

no 7 J=.. 

LA(I J, I, J)=LI( J) 

CONTIUUr 
L A A ( I R » J. f - o 

LAA( IRt . o/A+i’TL+M i/A-»*2 

L A A ( I R T ' T A ) = o 

laa{ IR » 2,r:) =-’( b: +; c )^'B3/A«-»z 
LA A{ IR t » 2 ) = ■ o 

LAA{ IR t'Ati )=-( 32+/- ) '^ {■» 1® + N /I 'io + ( N-Be )-»B3/3® ) / h*^2 
LAA(IR»4,.J= c 

LA A( IR 1 4 1 2 ) =*• ( 3 , +y.': ) •> ( “ Is +N / i j o +N*B3/3o ) / A -**2 
LAA»( IR 1 4 , 2 ) = 0 
LAA( IR,:'t2)= ® 

LAA( IR,l.3)= c 
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L4A( P c 

LAA( P , 3, : )=' o 
L A A ( J F ) '■< t — c 
LAA( )= c 

09 Pr: JI = SIK-j. 

IJ = IR+;- JT 
DC 9 I = , j 
S C “ : ' c 
SD= o 

3 =. f 

CO J = 

=L i { I J , J , ) •■■-L ! J , I , J 5 +S E 

SC=SC+L.-.A( IJ , J , )-'LaiJ,I,J) 

SD=s c+L aa I j , j , ^ } -L . { I j , T t J ) 

LAA( IJ- . ,I,?) = S -.''■|(!J,T5 
LAA( IJ‘ , I , 1 ) = S:'. 

9 LAA( IJ--, fl,2) = :i 

S ^ “i./ «j 

L4A{ IJ- : i ) = 

LAA( IJ‘ : J 2 )= c 

LAA( IJ* , 2 )=*:'nLT( I J- : ,2) 

LAS( I j 

LAS{ IJ,,.':)=SC 
LAS(1J,::1)=S2 
22 CQNTIMU'- 

SC=:"P 
SD=a 
SF="' « 

on 1 1 j = '■ j ^ 

s’;=SP-+L4ft (1 , j, a--aM ,",ji 

SC=SC + LAMi,Jr D^L' (i,9, j} 

11 so=SD+L'M(:s j, a.‘L.a S5» j) 

LAS(.lta=3C 

LAS( 3, 

LAS( i,:a=s.: 

0 3 1?. !J = ' , P 

ss=- « 

SC=' O 

0 “i o 

D3 13 I = ',3 

SF=Sc+FU( IJtD-^LJ.aiJ.It'a 
SC=SC+Fli( IJ,I )-a.4MI J,! ,1} 

13 SD = S0 + FU{ I J,I )^^L4A fl J,It2) 
nnijti.t ■ )=«'3C 

Did J, ..,3) =-'10 

Did j, '■ J ")='■ a Muaa-S'; 

Did Ja:,.; )=LmS{ij,: ) 

01 d J, , ) =LA.^ ( J J,.' ) 

12 Did J» Vt';- )=L4Sd J,' )*-:'JNd Jaa 
ICJ=IR 

IC2=l 

IC''=I 

65 DO l‘'i 1 = ' I-'': 

on 1 A j=: , ; 

14 l.)I (2 1 1 ! J ) = LS { ;. T I ad«-0I{ ..t5, J) +FUdtn*DI (1 16 , J) 
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0-1 ij = :,ir 
in c I = . ^ -1 
5C=- . 

$£='.: ■ : • ; 

!? v= -> ' ' ' f 

O'] j = , ,:f 
5F=L'{ T J, I, J j 
oc=sr + ii { ij, J, ) .vS"" 

'TO=GO+:;l ( ■■)*SF 

G •=:■ ■ + :'T { ij , J, n #SF 

x7 C-JlT ruj.. 

SF = FU( IJ, 'I ) 

Di{ij+ , ; , )=‘',c+ni ( I J,-'., : )-'^5F 

0I(IJ+,,i,l)=Gl+ll(lJ,!S,F)-5»'SF I 

DI(IJ+ 5 I , - + 01 ( !J,4,^ )^SF 

IF CaHTPUJ'-- 

5.4(1 1’ )=!)I{ IF+l ■.) 

ami,.' )=i)n 1?, + .,, ; 

AACi," ni(Tr. +,'.,1,13 
sc= ^ 

S0=.. c I’ 

se=D^ [ 

CO 16 i 

sc=sc+oi( tD^jDDd) I 

SD=SD+CI{ IP4: , I, Mi^CCCI ) 5 

S- = S-+D1 ( IP + \, !, M^CCd ) I 

16 CONTIMU'- 

AA (2 ,1 )=SC I 

AA )=sc 

AA(2, ■ )=■ G- 

CT=Ad:., ,)'dM i,23--AAd,2)*»AA(2a) 

I F(DTc w-.' , ) GO TO : 3 

NUd ) = {.^.4( '‘AA(..,,::)-AA(2,3)-»AAa,?))/CT ! 

■3U{?.)=('.AC.,' ) ‘AA( )' AA(2,l)*AA{l,5))/DT ' 

DO 2 IJ = „, IP, I 

niaiK ij,-.., 3->mjc )+ordj,5,:)»NU(2)+Didj,5,33 I 

nwc=01 ( I,J, A , )*nu( .d+DI dJ,6, ,:)«NU(2>+CI d J,6,3) I 

'riN(Tj)=DT. I 

W I K IJ ) — iJ W L ■ 

.;:A CniT PiU„: ! 

IF{ IC?,- '-Jo ) R-'.TUP.K , : 

L) 0 t-* 1 J ~ I r . i 

I F(TT{ i J) P.Sdi V-'. „) GO TO 61 

IF(TT( TJ )o ■ C. Q ) GO TC 62 I 

GO Til . ; 

cl . IFdI0(IJ).-,L'' o . ) GO TO £ . 

GO TO ^'1 : 

•..3 IF(TirM IJ )oG"'o ,:) GO TO 6:. i 

oP Di( IJ,",: )=>:.. ; 

Dl d J, )=n« ■■ I 

Did J, G,d=' 

TC" = jc:.' " ' ■ . ' : 

6 COOTIOtl,. . • , 

IFdd.d.'.cIR ) ddURO i: I 

TFdCP G :d GO TO 66 ; I 
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IC3=> 

RETURN 
c4 IC2 = ‘' 

GC T'3 f .. 

IS PRINT'S. 

2 FnRNAT(. T’T.viNtNT IS ZERO<f ) 

STOP 

rND 

IBFTC SUBS 

SURROLliri* C3J •{TT,TFUT,XINP,ABI,AB2,T> ,SC) 

C 

C THIS SUBft'ILT INS CsLCULATI: THE OBJECTIVE FUKCTIOM 

C 

COMMON /?,L.A/ 

CO^'MQ^j /BLSl/ .) 

COMMON /BLil/ ) ,Cl»C2tBl,B2»B3»DlV,lRtIR1 

COMMON /BL*'! / NIN 

niMON-SIUM TT ( 'H , WC { ' ) j KIN{ *=■ ,6 ) ,TaUT( 5) ,XlNRf5 ,4) 

AX=e 456 «' 24 o «34 c/(., ^ 

SD=AD(1,? j-^ITTd)^* •?' T ) ;2.= +AD(2tl )*(TT (I )‘ T‘ ) 

•Sc=So 

SC=WC( l)«l^!o/( 

CO 45 IJ=2,IR 
SC=SC+WC( IJ 

SD=S0+AD(i,l)»(TT(IJ)^Hi-,r;-TCLT{IJ-i)«»2)/2o+AD{2,l)*{TT(IJ)»T0UT( IJ 
D) 

D Q 4 6 I ~ j, 1 4 

46 SE = SE+( AD( 1, 1 + i )■»( "^T ( IJ )•»*'> TCUT ( IJ-,l )«*2) /?.*+AD{2, I+l )*{TT( I JI^-TQ 
1UT(I J- ■ ) ) )*)!IMR( IJjl ) 

45 COMTINUr 

sc=sr> 'S: Co 972/2''-3g.. 2 
sc=sc+(sc+sE 

SC =S C* *’ o . 

RETURN 

FND 

IBFTC SUB? 

SUBROUTINE CBJ' n{TT,TCUT,XINR,A3l,AB2,T< ) 

C 

C THIS SUBROITIN- IS TO CALCULATE THE DERIVATIVES OF THE OBJECTIVE FUNCTION 

C 

COMMON /Bin/ AD{2,:>) ' 

COMMON /8Lgy ,A3C 

COMMON /3L11/ DD(4),CI,C?,B1,B2,B3»DIV,IR,IR1 

COMMON /BL14/ WC,A 

COMMON /6L.2*"/ J.lUT 

COMMON /RLJi/ ''IN 

DIMtNSlU.'l XCUT( ';,6) 

OlMEILSinM IT ( S ) , i(^C { .S ) ,XIN { 5 »6 ) »TOUT( 5 ) tXINR { 5 ,4) 

on 41 I =.,.,A 

4,1 XIN( It I)=; . s 

AX=o^56#?.4o^:‘6> ./(IS ..'U^’o^ABa) 

XIN( :i,M={AC (1,1)«1T(1)+AD{ :,i) )*AX 
XIM( it 5) =XIN (I t'-I'^ABC 
X I N ( 1 1 ) =1 4 (, / { o 9 .A P 1 s.' AB 2 ) 

XIM( It U) =/IN (1 £o 872 / 2439«2 



